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Abstract 



A topological theory of d-wave superconductors is derived in this the- 
sis. Ginzburg-Landau theory describes superconductivity by defining a 
complex order parameter and applying Landau's theory for phase transi- 
tions. However, there is no local, gauge invariant order parameter for a 
superconductor and classical order is no appropriate description. Topo- 
logical order has proven to be a powerful tool for describing the Quantum 
Hall efi^ect and it gives also an appropriate description of s-wave supercon- 
ductors. For d-wave superconductors there are gapless excitations at four 
points on the Fermi surface. The topological theory for superconductors 
exhibits a similar structure for the ground state as found in the s-wave 
case. However, the gapless excitations destroy the topological degeneracy 
of the ground state and introduce an additional degeneracy in one of the 
flux sectors. In search for regularities which reflect the topological order, 
I include a magnetic point impurity and focus on the effects it has on 
the Casimir energy. Indeed, it is shown that the energy shift reflects the 
topology to some extend. However, this is work in progress and further 
computations will be done. 
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1 Introduction 



1.1 Classical and Topological Order 

In classical physics we can describe an N-body system by a probability distri- 
bution P(ri...rN). If we change the external conditions a little, the distribution 
will typically change continuously. However, if we change too much it might 
cause a discontinuous change in the properties of the system, a phase transi- 
tion. Ginzburg and Landau developed a very successful theory that describes 
phase transitions in terms of a (local) order parameter which is zero in one 
phase but finite in the other. Order is a property of a (broken) symmetry of the 
Hamiltonian. The ground state breaks one or more symmetries of the Hamil- 
tonian. Each broken symmetry generates a gapless particle, which is called 
Goldstone boson. There is, however, a loophole for this mechanism. In pres- 
ence of a gauge field, the Goldstone boson will combine with the gauge particle 
and form a massive, spin-1 particle. This is known as the Higgs Mechanism. 

With help of the order parameter 0, one can write an expansion of the 
Landau free energy in terms of (f>. This expansion will look the same (except for 
the coefficients) for systems with similar symmetries. Therefore, systems with 
similar symmetry properties have the same critical behavior and are said to be 
in the same universality class. 

Although this theory is very successful for classical phase transitionsQ it fails 
to describe the Quantum Hall Effect or charged superconductor^ In both 
cases no local, gauge invariant order parameter can be found [1]. Therefore, 
it is more suitable to describe these with topological orders to explain their 
characteristics. Topological order is a property of the ground state wave func- 
tion 7p{fi...riq). Classical order is a function of the probability distribution 
P{ri . . .tn) = |V'(^i • ■ • ^^f)P (at finite temperature). It misses the phase of 
the wave function, which is a quantum effect. Therefore, it can only partially 
describe the internal structure of the ground state. 

Topological order is robust against local perturbations as for example im- 
purities. The ground state degeneracy for example is not a consequence of a 
symmetry of the Hamiltonian. Therefore, impurities, which in general destroy 
the symmetry, do not have an effect on topological orders. However, it is sensi- 
tive to the topology of the underlying manifold. For trivial topology the ground 
state is unique. The simplest example for a non trivial topology is described 
by a hole in the manifold. A closed path around this hole cannot be deformed 
continuously into a point, it is said to be non-contractible. The paths can be 
labeled by how many times they wind around this hole. This will lead to degen- 
erate ground states. Thus, topological orders can be characterized by analyzing 
the ground state degeneracy [2 

Another probe of topological order is to study its defects: the low energy 

^However, even at the classical level symmetry breaking cannot describe all kinds of phase 
transitions, see for example the Kosterlitz-Thouless phase transition 

■^Although Ginzburg-Landau theory, in principle, does not qualify for superconductors, it 
gives nevertheless an accurate description of for example the H-T phase diagram 
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excitations, so called quasiparticles. These can carry fractional charge and have 
fractional statistics. In contrast to the Goldstone Bosons of the classical order 
these particles are normally massive. Fractional statistics of point particles is 
generally a feature of two dimensions. Consider two identical particles. Taking 
one particle around the other, the wave function will acquire a phase. In three 
dimensions this path can be continuously deformed into a point. Thus the wave 
function cannot be changed by that path and the phase is a multiple of 27r. 
Therefore one can only obtain two different values of the phase for exchanging 
the particles. For bosons this phase is zero corresponding to integer spin and 
Bose-Einstein statistics. For Fermions this phase is tt with half-integer spin and 
Fermi-statistics. In two dimensions this argument is not valid anymore. The 
encircled particle prohibits the deformation of the path to a point. Thus the 
phase can in general be arbitrary. The particles are therefore called anyons. 
They interpolate continuously between fermions and bosons. 

As the wave functions are not single valued in two dimensions, one can 
count how many times one particle circles around another. If the underlying 
manifold is not simply connected, there are also closed path on this manifold 
which cannot be continuously deformed into a point and will put a phase to the 
wave function. This forms the so called the braid group. For two particles on 
a simply connected manifold it is just the group of integers. If both particles 
are identical exchanges are also allowed. This is described by including half- 
integers. For N particles on a non-simply connected manifold the braid group 
is generally more complicated and does not have to be abelian [2] . 

If one has a finite system with a boundary, there will also be edge excita- 
tions. Their structures are extremely rich and, in general, they provide a more 
complete and a more practical measurement of topological order than the bulk 
excitations. To get a better idea of these excitations consider a finite droplet 
of a FQH liquid. Although it cannot be compressed changing its shape will not 
cost much energy. Thus one can visualize edge states as the surface waves on 
these droplets. In this case the edge states will be gapless, [5], in general they 
can also be gapped. 

The concept of topological order was first introduced in a study of spin liq- 
uids. The first experimentally observed state with non trivial topological orders 
was the superconducting state. In standard literature superconductors are of- 
ten described by breaking of a gauge symmetry. The gauge photons will acquire 
a mass. Thus, the magnetic field will be excluded from the superconducting 
material. The order parameter is taken to be proportional to the pair wave 
function {ca{k)c-a{—k)) and is in general a complex- valued function. However, 
it is not gauge invariant and, therefore, not a suitable candidate for the order 
parameter. On the other hand, gauge invariant objects which could serve as 
an order paramter are not local. So there is no local order parameter in the 
superconducting state [T]. The main success of topological orders came with 
the FQH liquids. Although not all of the FQH states can be explained by it, it 
provides a good description for simple so called filling fractions v = —. 



7 



1.2 Outline 



This diploma work analyzes whether a high-temperature superconductor can 
be described by a topological theory. It has already been shown that this is 
an appropriate description for conventional superconductors [T] . The supercon- 
ducting state exhibits various features of topological order, i.e. its excitations 
are gapped and it is sensitive to the global topology of the underlying mani- 
fold. In high-temperature superconductors, the gap function is not of s-wave 
but of d-wave symmetry: Ak ^ + ky. In contrast to the s-wave symmetry 
of conventional superconductors, it vanishes at four lines in momentum space. 
At the points where these lines intersect with the Fermi surface, the energy de- 
pends linearly on the momentum and the excitations are gapless. It is an open 
question, how these gapless modes affect topological properties. 

In the first part of this thesis, I derive the low-energy theory, starting from 
the BCS Hamiltonian. It is shown that the appropriate description is a Dirac 
Lagrangian for free, massless fermions, coupled to a topological gauge potential. 
This model is then probed by taking a torus to be the underlying manifold and 
the results are then compared to the ones obtained in the s-wave case. Although 
the quasiparticle fields couple to a topological gauge potential, the model itself 
yields results, which differ in some important ways from conventional topological 
models. 

Some aspects of topological order are destroyed by the presence of gapless 
modes. It can be shown, that the energy spectrum becomes sensitive to local 
perturbations. The question arises, whether these modes destroy all topological 
features or whether one can formulate some weaker statements, which also hold 
in the presence of gapless quasiparticles. More specifically, in the second part I 
look for other features which depend solely on the topology. It is argued that the 
energy shift under local perturbations may be such a one. Although the ground 
state energies arc not constant under these perturbations, there is evidence that 
they at least behave in the same way. It can be shown that the energy shift 
vanishes for all ground states in first order perturbation theory. The expressions 
for the second order contribution are derived and regularized. For symmetry 
reasons, it can be argued that some of them show similar behaviour even in 
second order. 

In the introductory section, I will give some necessary background material. 
I begin with a short review on topological field theory. In particular, I give a 
short introduction in BF-theory and state two applications, the quantum Hall 
effect and superconductors. I proceed with a section on the Casimir energy and 
give some comments on regularization in quantum field theory (QFT). 

1.3 Topological Field Theory 

A field theory is called topological if neither the action nor the observables 
depend on the metric of the underlying manifold. The easiest way to achieve 
this is to require both quantities to be explicitly metric independent. Examples 
are BF-Theory and Chern-Simons theory, both of which can be solved exactly for 
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simple cases. The observables for Chern-Simons gauge fields A are Wilson loops 

e" f ^ . In this section, I start with describing BF-Theory and its application, 
first in the Quantum Hall EflFcct and then in the superconducting case. There 
is also a discussion, how the latter behaves on a torus. 

1.3.1 General Chern Simons Theory 

In contrast to BF-theory which can be formulated in any number of space time 
dimensions, Chern-Simons theory (CS) only works in an odd number. Mostly 
one is interested in d = 2 + 1, where the CS Lagrangian is quadratic in the fields. 

C = ^e^'^a^d.ax-a^r (1) 

where k is a positive integer. 

This Lagrangian is gauge invariant, if one can neglect boundary terms. This 
is for example possible on a compact manifold such as a torus. Varying ^ with 
respect to a^, one obtains the Euler Lagrange equation. 

Ke^'^^F^x = r (2) 

In contrast to Maxwell's equations there are no derivatives involved. The cur- 
rent conservation is ensured by the e'^''^ tensor. In particular the lack of time 
derivatives implies that the fields are non-dynamical [6]. Without a current, the 
only degrees of freedom are topological ones. Also the different components of 

do not commute with each other but form a canonically conjugate pair. 

To see this, first note that ao can be regarded as a Lagrange multiplier as 
there are no derivatives acting on it. Thus it constrains Ui to fulfill: 

e^'d.a, - f (3) 

Explicit solution will depend on both, the topology of the manifold and the 
density j^. Calculating the Euler-Lagrange equations for the remaining terms 
one finds that 02 is the conjugate momentum to ai, leading to the following 
commutation relations: 

[ai(f),a2(y)] = -S{x-y) (4) 

K 

As Qj is a gauge field and thus not gauge invariant, it cannot serve as an 
observable. These are given by the Wilson loops around non-contractible paths 
7. These are gauge invariant, global operators. 




1.3.2 Quantum Hall Effect 

The best known example of a topological ordered physical system is the Quan- 
tum Hall eflfect (QHE). Both the Integer Quantum Hall Eflfect (IQHE), discov- 
ered 1980, and the Fractional Quantum Hall Effect (FQHE), discovered two 
years later, had a great impact on condensed matter physics. 
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Classically the Hall resistivity should be a linear function of the electron 
density. In the quantum regime, i.e. at low temperatures and high magnetic 
fields, the resistivity is not linear any more but develops plateaus at certain 
"magical" values of the filling fraction v, which is the corresponding quantity 
to the classical density: 

Number of electrons /?$ 
Number of states N B 

The IQHE describes the plateaus for v being an integer number, fractional fillig 

factors can be described by the FQHE. Moreover, states with simple rational 
filling factors, i.e. v — are more stable and easier to observe then states with 
more complex ones, i.e. v = ^. These require lower temperatures and cleaner 
samples. The IQHE on the other side is visible even for rather dirty samples. 
The observables are also insensitive to the microscopic details of the sample. 
By measuring the Hall resistivity one obtains a very precise value for the fine 
structure constant a. If the resistivity depended on e.g. the physical dimensions 
of the sample, a measurement of that accuracy would be impossible. But the 
quantization is not the only interesting property of this new state. Because of 
the light electron mass, the 2-dim electron gas cannot condense to a crystal but 
will instead form a quantum liquid. However, it is much more rigid than a solid 
as it can not be compressed. Furthermore, though it does not develop a spatial 
order, the movements of the electrons arc highly organized. 

I want to consider a two dimensional electron gas with an additional magnetic 
field perpendicular to the plane. Classically one would say that the electrons 
move in circles. In quantum mechanics the energies of the electrons will be 
quantized in terms of the cyclotron frequency Wc, so that each electrons can 
have an energy of E = {2n+l)JkVc- These levels are called Landau Levels (LL). 
The Lowest Landau Level (LLL) contains all the states with n = 0. Taking the 
spin into account, it is filled at = 2. With increasing magnetic field both the 
number of states and the energy difference of the Landau Levels will increase. 
So for high enough magnetic fields and low enough temperatures one can always 
restrict oneself to the LLL and ignore the higher energy levels. 

As stated earlier the QHE is described as a quantum liquid. In the IQHE 
all the states of the LL are filled. If one tries to compress the liquid one would 
have to excite electrons to a higher energy level. However this requires a finite 
energy due to the energy gap between the LL. This will prevent the compression 
and thus lead to the very low compressibility. 

For the FQHE at filling fraction = ^ = 2p+i explanation is a little 
bit more subtle as, in that case, there are vacant states in the LLL. The energy 
gap is in this case provided by the large Coulomb repulsion of the electrons. 
The system can be described as neutral plasma which will try to spread out the 
electrons as homogeneously as possible. By compressing the plasma one has to 
overcome this large Coulomb repulsion. 

However the system can also be described with help of composite parti- 
cles, i.e. particles to which flux quanta are bound. This can be achieved by a 
Chern-Simons field. If it binds 2p flux quanta to each fermion, it will become a 
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composite fermion, moving in a reduced effective magnetic field: 

B = B- 2p<i>p 

Thereby it will fill the LLL and the effect can be described as the IQHE of 
these composite fermions. If it binds 2p + 1 flux quanta to the electron then 
a composite boson will be formed. These bosons will experience no effective 
magnetic field at all and therefore will condense into a Bose condensate. This 
condensate is then incompressible by a many-body Coulomb effect. |5] 

In addition it's interesting to see what happens when we insert or take away 
an electron: the system at filling = ^ will react as if m quasielectrons resp. 
quasiholes have been inserted, each having a fractional charge of e* = and 
also fractional spin. This also emphasizes again the fact that the appropriate 
description for the QHE is by a Chern-Simons theory. 

1.3.3 BF-theory for superconductors 

In my discussion I will use a Chern-Simons like theory, usually referred to as 
a BF theory. As I also consider vortices, two Chern-Simons fields are needed. 
The name is simply reflecting that the second gauge field normally is denoted 
by b. The two fields will couple to two different currents. In my case will 
couple to the spin current j^, while will couple to the vortex current j^: 

C = -e^^'^d^ax-a^f -b^^r (5) 

TT 

The coefficients in the Lagrangian are chosen to give the correct braiding statis- 
tics. When a quasiparticle goes around a vertex, the wave function should pick 
up a minus sign. In this case the different components of commute with 
each other, the same for 6^. But one obtains nontrivial canonical commutation 
relations between Oj and bj: 

[ai(f), 62(27)] = iTT5{x-y) 

[02 (y), 61 (2/)] = -iT:5{x-y) 

Calculating the Euler-Lagrange equations for the Lagrangian gives the cur- 
rent conservation for both currents: 

IT 

r - -e^^^d.ax (6) 

TT 

After separating the imaginary time and space components and doing some 
partial integrations, the BF action can be written in the following form, 

S = — j d'^x e^^ dibj + aoe''^ dibj + boeijdiQj (7) 
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Analogously to the discussion of the general Chern-Simons theory, oq and 
60 can be regarded as Lagrange multipliers and the solutions for and bi have 
to fulfill: 

TT 

f = -e^'d^a, (8) 

TT 

1.3.4 BF theory on the torus 

By imposing periodic boundary conditions, one turns the Euclidean plane into 
a compact manifold, a torus. It can visualized by taking a parallelogram and 
identifying opposite sides. It's straightforward to solve ([8]) on a torus: 



b^ = d,Ki + ^ (9) 

where Ag and At, are periodic functions on the torus and on and bi are spatially 
constant. After reinserting ([9]) into ([7]), the action is given by: 

S }l [ d^^e'^{dt{d,Ka)d,Kb + ^d,Kb + d,ka^ + ^) (10) 
TT J Li Lj LiLj 

The second and the third term vanish, as they are total space derivatives of 
periodic functions. The first term can be integrated by parts and written as 
e^^ {dtAa){didjAi,). This term vanishes also, as the derivatives commute with 
each other. The final Lagrangian is given in terms of the spatially constant 
functions a,j and 5^: 

TT LiLj 

The quasiparticle and vortex number are quantized. This requires the corre- 
sponding gauge fields to be compact: di and bi to be angular variables defined 
only modulo 27r. In a superconductor, the flux is quantized in multiples of tt. 
Therefore, both gauge field can have only two different values, (for no flux 
quanta present) and tt (for one flux quanta present). 

One can also approach this problem with algebraic methods. For this pur- 
pose one defines operators which generate non-contractible loops. With only 
one type of quasiparticle present, there are three such operators, two which 
take the quasiparticle j around the two cycles of the torus, Tj and pj, and one 
taking quasiparticle j around another one, aj. Any process can be written as a 
product of these basic ones, aj commutes with the other operators, but Tj and 
pj don't. Assume you want to measure the flux through both cycles. By taking 
the quasiparticle around one cycle, you measure the enclosed flux. However, 
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at the same time you create a flux through the other one. Therefore, the two 
measuring processes don't commute with each other: a result of the non-zero 
commutation relations between the different spatial components of the gauge 
field. 

With two types of particles the braid group becomes more complex. In 
particular, one has much more operators to deal with. Therefore, I restrict the 
discussion to just one quasiparticlc on the torus. Then the only non-contractible 
paths left are the ones around the cycles of the torus. 




Figure 1 : A vortex tunnc;ling process inserting a unit of magnetic flux inside the 
torus. In this visualization it also leaves a flux loop outside. But as electrons on 
the surface on the torus cannot wind around this external flux, we can neglect it 
in the further discussion. This tunneling process connects ground states labelled 

by opposite values of the Wilson loop e^c "^'^ = g**^*** where $m is the magnetic 
flux threading C. 

For each gauge field you can define two corresponding Wilson loops given 
by: Aj = e^"'^ and Bj = e^^^ . Acting, for example, with Ai can either be seen 
as measuring the a-flux in the j-direction or as inserting a b-flux in the other 
direction. Each set of eigenvalues for the Wilson loops A^ and Ay defines a 
sector. They can be connected by tunneling processes: a vortex antivortex pair 
is created and recombines after having moved along eg. the non-contractible 
loop in x-direction. This process creates a flux quanta in the y-direction and 
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moves the system from one sector to another. 

In the fully gapped case in an s-wave superconductor both types of quasipar- 
ticles are massive. Thus, tunneling processes are suppressed and the different 
sectors are separated by an energy barrier. Each sector has a unique ground 
state, so there is a four-fold degeneracy of the ground state on the torus [T]. 

1.4 Casimir effect 

The Casimir effect is a manifestation of the zero point fluctuations of the vac- 
uum. The simplest example of is two neutral, conducting plates hanging parallell 
to each other at a small distance. Classically, there is no interaction between 
the plates, but quantum-mechanically, the electromagnetic fields oscillate even 
at zero temperature and cause an effective, small, but measurable, attractive 
force between them. The size of this "Casimir force" was first calculated by H. 
Casimir. The total zero point energy (per unit volume) of the electromagnetic 
field is infinite, simply because there is one harmonic oscillator for each point 
in space and their zero point energies are all summed up. Nevertheless, the 
difference in zero point energy for different positions of the plates is finitelfft is 
this finite energy difference that is manifested in the Casimir force. 

In a more general setting, the Casimir force can be thought of as the finite 
response of an (interacting) field theory to varying boundary conditions. Also, 
there might be couplings to external potentials. In our case, 2d fermions are 
coupled to a (topological) gauge field in the presence of a point like impurity 
potential. Since the theory is defined on a torus (Li, L2), the relevant quantity 
is the energy density, £{Li,L2)- 

In order to calculate the Casimir energy, the ultraviolet (UV) divergences 
in the sum of the energy levels in the Dirac sea, (or in the case of bosons, the 
sum over the zero point fluctuations) must be regularized, and the interactions 
renormalized, in order to extract a finite and physically meaningful result. The 
actual procedure will be briefly outlined in the next subsection, and described 
in detail in section 4.4 below. Since the full description is rather technical, 
however, I will first give a non-technical summary of the logic: 

1. Since the gauge field is topological, it does not have any local degrees of 
freedom. More precisely, there are no fluctuations in the topological gauge 
field which could contribute to the Casimir energy. 

2. I will neglect the coupling between the fermions and the gauge field, except 
for the coupling to the external fiuxes threading the torus. Technically this 
means that I will only consider diagrams with a single fermion loop. 

3. I will assume that the impurity potential is weak, and can be treated in 
perturbation theory. The leading order term thus correspond to a free 
fermi theory coupled to external gauge fiuxes defining the various sectors 
described in the previous sections. 

•^This is true for an electromagnetic field in the presence of perfect conductors, but in a 
more general context, the Casimir energy does not have to be finite. 
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4. There are several ways to calculate the Casimir energy of the nonintcract- 
ing theory, but I will use a method that will generalize to the interacting 
case. I relate £ to the full interacting fermion Greens function G{x, y, u) 
on the torus, and the free field case amounts simply to use the free Greens 

fimction Go{x,y,u)). 

5. In order to isolate and subtract the leading free field theory divergence, I 
write the full torus propagator as G{x, y,ui) = Go{x, y, w) +^^_^q Go{x + 

rh,y,u}), where the sum is over windings (in multiples of Z/(27r)) around 

the cycles of the torus. This representation isolates the short distance 
singularity in the first term (corresponding to the free case) which can be 
directly subtracted as described in section 4.2. 

6. Interactions with the impurity potential introduces new singularities, and 
in order to cancel these, the impurity interaction has to be renormalized 
- the procedure is described in section 4.5 . 

7. The respresentation of the full propagator as a sum over windings is valid 
only when the sum actually converges. As is shown in section 4.5 this 
is not the case for w = in one of the flux sectors. This introduces a 
spurious infrared divergence, and the associated technicalities are explaind 
in section 4.5.1 

1.5 Regularization and renormalization of the impurity 
potential 

It is a well-known problem in QFT that point like interactions lead to divergent 
integrals. Those have to be renormalized in order to extract physical quantities. 
In the second part of this work, I consider a magnetic point impurity, which is a 
singular potential and causes similar problems. The quasiparticles scatter at the 
impurity, which shifts the energy levels. We are interested in the sum of all these 
shifts. However, as discussed above, this sum is in general infinite. To define the 
expressions we must first regularize and then renormalize the interaction. The 
divergencies first arise in second order perturbation theory, where the particles 
scatter two times at the impurity. If one considers this process in free space, it 
is well-known that it gives a divergent contribution which must be isolated and 
subtracted in order to calculate the renormalized Casimir energy. The cause 
of the divergence is the singular nature of the impurity. It scatters into all 
momentum states with the same strength. When the particle scatters at the 
same point, the potential does not provide a natural momentum cutoff for the 
integral over the internal momentum. This is nothing but a loop in QFT and 
we can use the same techniques to regularize it. As soon as the particle travels 
along a non-contractible path between the scattering events, the minimal length 
of this path provides a cut-off and makes the expressions convergent. 

There are various ways to provide a momentum cutoff A. The detailed 
form of the regularization does not effect physical results as long as it does 
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not break any symmetry of the theory. Here, the regularization is achieved by 
point-spHtting, i.e. by considering the interaction to happen at two different 
points afo and xq + 5x. As long as 5x is considered finite, though very small, 
the expressions are finite. Taking the limit, 5x — > 0, restores the original theory 
and gives rise to divergencies. In order to cancel the divergencies, I add local 
counterterms in the Lagrangian. As we are interested in topological features of 
the model, it is important that we don't have to introduce nonlocal counterterms 
in order to regularize the expressions. 

Consider first the case where the quasiparticles are massive. In order to 
distinguish between the different topological sectors, the quasiparticles have to 
travel along non-contractible paths. On a torus, these are the paths around 
the two handles. The different topological sectors are given by the flux quanta 
through the holes of the torus. When the particles are massiv^l, each path 
around a hole is suppressed by a factor e"™"^, where m is the mass of the 
particle, L = {Lx,Ly) the torus size and n the winding number. Therefore, 
windings are suppressed in the limit of a large torus. As a result, the ground 
states in the different sectors are degenerate. If a local impurity is inserted, the 
quasiparticles scatter on this impurity. However, as long as the windings are 
suppressed, a local impurity can not change the degeneracy of the ground states 
or, more generally, any topological feature. However, by inserting a non-local 
perturbation, windings could be facilitated and topological features could be 
destroyed. The same is true for nonlocal counterterms. 

In presence of gapless excitations, the situation becomes more involved. The 
exponential factor vanishes for zero mass and windings around the handles are 
no longer suppressed. As a consequence, the ground state degeneracy is lifted. 
In contrast to the case where all quasiparticles are massive, even a local impurity 
can induce different energy shifts for the ground states in the different sectors. 
The issue addressed here is whether the energy shift shows some regularity. 
Also in the case of gapless excitations it is important to just consider local 
perturbations and to regularize them with local counterterms. 



*In this work, I use both the notion of mass from particle physics and the notion of gapped 
particles from condensed matter interchangeably. 
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2 Dirac Theory 



2.1 Quantization of Dirac Theory 

In the superconducting phase the electrons form Cooper pairs. To break these 
Cooper pairs, energy is needed. In normal superconductors this energy does 
not depend on the direction in fc-space: the gap function has s-wave symmetry. 
This is different in high-Tc superconductors. There the gap function has d-wave 
symmetry 

A{k)=VAikl-kl) 

So along two Hues in k space, k^ = ±ky the gap function will be identically 
zero. Thus there will be four points on the Fermi surface where the excitations 
are massless. I will denote these with fcj , i = 1, ..,4. 

ki = ^{1,1) = -ks and ^2 = -^(-1, 1) = -^4 







V 







Figure 2: fermi surface in d-wave superconductors 

First the quasiparticle fields are expanded around around the nodes to ob- 
tain the low-energy description. Because of the d-wave symmetry of the gap 
function, the obtained model will be linear in both the time and the spatial 
derivatives. This indicates that Dirac theory is the natural language to describe 
these excitations. The quasiparticle fields can be seen as massless fermions, 
which are coupled to a gauge field A^. Therefore I will give a short review of 
two-dimensional massless Dirac theory. I start with the massive case and then 
discuss the massless case in detail. 

Normally two dimensional 7-matrices are sufficient to describe this theory. 
The most convenient choice is to use the three Pauli matrices as Dirac matrices. 
However, I will need to use four-dimensional matrices. 
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After the expansion, there are two independent fields at each node, for spin- 
up and spin-down respectively. If the particles at every node act independently 
of the other nodes, then a two-dimensional description would be possible, con- 
sisting only of spin-up and spin-down particle at a given node. However, the 
nature of the interaction is to couple particles with opposite momentum and 
opposite spin. Thus, the spin-up particle at node 1 can not be treated inde- 
pendently of the spin-down particle at node 3. Due to this interaction, you 
need not consider an eight-dimensional analysis. The quasiparticles are created 
by destroying Cooper pairs. In the Hamiltonian there are two types of terms: 
cl{q}Ca{q) and c'^\q)c^^l^{—q). It is easy to see, that e.g. quasiparticles at 
node 1 don't interact with those at node 2 and node 4. So, the theory can be 
restricted to one pair of nodes. 

In BCS theory of superconductors, the conserved quantity is spin, not charge. 
The goal, therefore, is to find a formulation, where the spin current can be 
identified with the conserved Dirac current. The naive choice for the spinors is 
to put all fermion fields into the spinor and their adjoint into ^. However, it 
is quite easy to see, that this is not an appropriate choice. First of all, because 
it does not even enable you to write it as an Dirac equation. But if one wants to 
interpret spin as the conserved charge, the spinor needs to be an eigenstate of the 
spin operator. Therefore, you combine spin-up operators with the adjoint spin- 
down operators in the spinor, see [3]. Although 1 use a four-dimensional matrix 
representation, the theory itself is two-dimensional. Therefore, only three of the 
four Dirac matrices are needed. The third, 73 is redundant. It is still useful, 
though, to define it. It will appear later in the discussion as one of generators 
of the chiral symmetry of the two-dimensional system. 

Throughout the whole discussion I use the metric gf^" = diag{l,—l,—l). 
The Lagrangian for the spinor fields can be written in following form: 



C{x) = *((x)7''(z9^- A^)-m)*(x) (12) 

using natural units where c — h — 1. It is helpful to expand the spinor fields in 
eigenfunctions of the Dirac equation: 

^ix) - y'^^=^aX(p>-''-+5;V(p)e^''- (13) 

a* and b'^ are operator valued coefficients, and v'^ are the eigenfunctions of 
the Dirac equation for positive and negative energy and E = pq is chosen to 
be positive. The factor ,^ was introduced as the spinors are normalized to 

•y/2E^ instead of 1. The Hamiltonian is diagonal in the operators and b'^: 

^ = J ^T.^p(<K + K^K) (14) 

^ ^ s 

a^^ creates a particle with momentum P , 6''^ destroys one with momentum —p. 
This can also be interpreted as creating a hole in the Dirac sea with positive 
momentum. 
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2.2 Massless Fermions 



In the absence of a mass term, an additional symmetry arises, chirality. In 
3+1 space dimension, one can define a matrix, 75, which is self-adjoint and 
anticommutes with the Dirac matrices: 

75 = no7i7273 (15) 

For every spinor 'f, which is a solution of the Dirac equation, 754* is also a 
solution. As 7I = 1, its eigenvalues are ±1. Thus, the eigenstates of the 
Lagrangian can be chosen to have a definite chirality, they are called right- 
handed for chirality = 1, left-handed for chirality —1. 

In 2+1 dimensions, there are only three Dirac matrices instead of four. Anal- 
ogous to the three dimensional case, one can define a 73 by claiming that it has 
to fulfill the Clifford algebra. So there are two matrices, 73 and 75, which 
anticommutc with the Dirac matrices and one, which commutes, 7375. These 
matrices don't commute with each other. Instead they form an SU(2) algebra 
like a spin. So the chirality in 2+1 dimensions has the same structure as the 
spin symmetry. Each matrix can be linked to a direction and only one can be 
chosen to be diagonal. 

2.3 The Effective Lagrangian 

First I will give a short review on BCS theory. I start with a rather general 
ansatz for the Hamiltonian: 

Ti. = Hq + Hint 

q,(7 

+ ^ Vk{q,^)cl{q + k)cl^{-q + k)c-,{~^ + k)ca{ci + k) (16) 

Ca and cl are the usual fermion operators and obey canonical anti-commutation 
relations: {a[^{q) , a p{(f )^ = Sa.pSq.q'- The interaction term, Tiint is invariant 
under global charge U(l) and spin SU(2) symmetries and conserves both, mo- 
mentum and spin. In BCS theory two electrons with opposite momentum and 
opposite spin form a Cooper pair. You expect, therefore, a non-zero expecta- 
tion value for terms like the pair operator: Ccr{q)c-a{—q)- You also assume, that 
the fluctuations in the pair operator are small, so you are allowed to make a 
mean-field approximation by introducing the superconducting order parameter 
(or gapfunction) , A: 

^£(9) = l]t^;t(a,5')(ca(a' + fc)c-a(-g' + fc)) (17) 
q' 

You can now identify k with the center of mass momentum of a Cooper pair, 
and also define a spatially varying gap function by a Fourier transformation to 
real space: 
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k 

If the electrons form a pair with zero total momentum, the gap function is 
constant: A^(g) = A(q). With magnetic flux tubes present, Cooper pairs have 
non-zero momentum and the order parameter becomes a function of their center 
or mass coordinates, R. 

The Hamiltonian for our system is given by: 

+ |ct(q)A^~(g)cL.(-g) + /..c. (19) 

where the dependence on the center of mass is coded in A^. Furthermore, I will 
take the gap function to have d-wave symmetry. As I am interested in the low 
energy physics, I will expand the creation and annihilation operators in terms 
of continuous fields '0(g) for q values near the nodes and integrate out the high 
energy degrees of freedom: 



As we have restricted the momenta to regions near the nodes, it is legitimate 
to linearize both the energy and the gapfunction. 



£.{kj + q) - ^ = 2kj ■ q + 0{q^) = vpqj 

A(fcj +q) = A(i?)[(fcj- ^q^)2 - (kj^yqy f] 

= 2A{R)[kj,xqx - kj.yqy] (21) 

where qj is the momentum parallel to kj, qj — kj ■ q, and vf denotes the Fermi 
velocity. Inserting this relation into (fT6|) gives: 

^ = y E V'It (9)^ (^.- + 910,,T (91 + i'limikj - q)4'U (q) 



3 = i 



-V'].i(9)A^(fc, - g)0]+2,^(q) + h.c. (22) 
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2 

(q)+h.C. 

One can simplify this expression by introducing two 4-spinors 

*l = (4T'^J-+2a,V^]+2,T>ja) (23) 

and writing it as a matrix equation. Furthermore I wiU write the gap function 
as A(i?) = j^e^'^. For simphcity I assume that the only space dependence lies 
in the phase iip, since amplitude fluctuations are costly in energy. 



/ 1 



*I(^)vAg2 



( 



^e-'-'P 



/ 



^e-'-'P 



V 



1 / 



e^^ 

\ 

e'^ 
e-'-f j 



^2(9) (24) 



Now do a Fourier transformation back to real space. 



V',.i(r) = j (02e-^%,,(g) (25) 

Note the "-" sign in the second line of ([^5]) . This convention ensures that the 
spinors Fourier transform in a simple way: 

« = / (02 e''"*:' (26) 

We would have obtained the same low energy fields, if we had expanded the 
fermion fields in r-space: 



(27) 
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In general, the phase of the order parameter will be space dependent. In the 
Dirac theory the matrices are constants. So in order to rewrite p2|) as a Dirac 
equation, I will redefine the spinors and, thus, remove the phase from the matri- 
ces. If the phase 4> is well-defined at any point on the torus, there is no problem 
since it can then be removed by a regular gauge transformation. The phase is 
divided equally between the spin-up and spin-down fields and the transformed 
spinor fields are still single- valued. However, cj) is an angular variable, so gener- 
ally it is not single valued. Consider a magnetic flux through one of the holes in 
the torus. The magnetic field on the torus will be zero, but the electromagnetic 
gauge potential is not. Instead, it will wind around the magnetic flux. By going 
around the hole including the flux tube, the phase does not return to its original 
value. As the magnetic flux in superconductors is quantized, the phase picks 
up an additional 27r. If now the phase is divided equally between the two fields, 
they are not single-valued anymore. Thus, the naive choice how to transform 
the spinors does not work. I will denote that part of the phase as the singular 
part. 

There are different approaches to deal with that problem. In [4j the au- 
thor considers vortices in a superconductor. Those vortices are divided into two 
groups, labeled A and B. The spin-up fields are transformed with the singular 
part of the phase, which comes from the vortices in group A, spin-down with 
those in group B. Finally, you average over all possibilities to distribute a given 
number of vortices into the two groups to get rid of the anisotropy you intro- 
duced. In this paper, a similar but much simpler method is applied. Here for 
simplicity, no vortices are present. The singular behaviour of the phase arises 
only of the fluxes in the holes of the torus. In general, the phase has both a 
regular and a singular piece. The regular part of the phase will be divided sym- 
metrically to both fields, the singular part, however, only to spin-up particles. 
Thereby, I also create an anisotropy between spin-up and spin-down particles. 
This will become manifest at the analysis of the ground state degeneracy. I 
include also an electromagnetic gauge potential to include external magnetic 
fields in the analysis. The transformation of the spinor fields is then given by: 



(28) 



where 





(29) 



this is equivalent to transforming the spinors by a unitary matrix U: 



U{r) 



(30) 
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As a result of performing this transformation there will appear additional terms 
in the effective Hamiltonian. As in the Chern-Simons theory you can either 
choose to have multi-valued wave functions or you can keep them single-valued 
by introducing a Chern-Simons gauge field a^. I introduce two matrices: 



ai 

They correspond to the a matrices in the conventional Dirac theory. The addi- 
tional terms which appear in (j24p after doing the gauge transformation, can be 
rearranged into two current terms coupling to different gauge fields. The Chern 
Simons gauge field: 

= -\^djips (32) 

couples to the conserved Noether current, which in our model is not the em 
current but the spin current. The total phase, that is both the Chern-Simons 
field and the electromagnetic gauge field, couples to the cm current. The spin 
current is denoted by 2% while the electromagnetic current is Jem- ^ more 
detailed derivation of the spin current, is given in Appendix A. 

jy = -5'|t'Aa2'I'i - 'I'2^'Fai*2 
Jem,j = ^]VF (33) 

With this, we can define the covariant derivative Dj — dj — iaj and write the 
Hamiltonian in Dirac form, 

n = -^\{ivFDxai+ivADya2)'^i-'^liivADxa2 + ivFDyai)^2 

+ \{d,^ + 2A,)3i^ (34) 

The Hamiltonian and the spin current depend on two velocities, one being just 
the normal Fermi velocity, vf- The other, vai characterizes excitations parallel 
to the Fermi surface. In cuprates, these velocities are not equal to each other. 
This anisotropy has consequences on the symmetries of the system. A more 
detailed analysis can be found in the discussion of the symmetries of the model. 



^3 











CTl 
CTl 



(31) 
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3 The Effective Lagrangian for a d-wave super- 
conductor 



To analyze our model, it is convenient to use a Lagrangian description, where 
the time dependence is included. For later convenience I will also transform the 
effective Lagrangian to Euclidean space. I will also include a dynamical term for 
the A field, —jF^,^, and introduce a vortex current, parameterized by a gauge 
potential 6^: 

r = -e^^^d^ipax (35) 

TT 

I account for the phase fluctuations by also introducing a mass term: 

The mass \l will depend on the size of the torus and go to zero in the limit 
of |L| ^ GO I will also introduce 70. As we are in Euclidean space, the gamma 
matrices have to fulfill the Clifford algebra for Euclidean space: 

{7^,7.} = 2(5^,, (36) 
which leads to the following definition of 7^ 



70 



I 

1 



71 = 70 • 71 

72 = 70 • 72 (37) 
The effective low energy Lagrangian is then given by: 

= -A + ^iA^-a^f-h^r+'-e^^Xd^ax (38) 

4 Z TT 

+ *7^(a^-m^)*-(A^ + ia^,-a^)J,^„ 
I introduced the abbreviation: 

*7^a^* = +*iWir7l9l +VA 7292*1 + *2t^F7l92 + UA729l*2 (39) 

The conserved Noether current in our model is not the electromagnetic current 
but the spin current given by: 

r = *7^* 

(40) 
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Which in detail reads: 

f = *it;A72*l + *2fF7l*2 

(41) 

For a more detailed discussion of the spin current see Appendix A. 
3.1 Derivation of the effective BF action 

First I note that the regular phase can be absorbed by a regular gauge trans- 
formation and then I identify + = Sa^ where 5a ^ is a "small" quantity. 
In the obtained equation 



^ = ^Fl + ^Sa^,-b^r + -e^''%d.{6ax-A^) (42) 

4 Z IT 

the Sa^ field can be integrated out. In the final equation we identify again the 
field as the singular piece of the gauge field and set = a^. 

4 TT 

+ *7^(ia^ - a^)* - hl{^e'^-%K + - JLf 

Zi TT 

1 \ 2 

+ i^e^''\d^b.){j^-Jem,^.)-^0^.-Jem,^f (43) 
TT ^ 

In order to simplify the analysis, I will assume that the vortex current van- 
ishes: = and I will keep only lowest order terms 



£ = ^e^^^fe^S^a^ + M/7^Z)^M/ (44) 
3.2 BF-action on a torus 

First I will discuss the BF part in Minkowski space. It is illuminating to rear- 
range the terms in the BF-Lagrangian: 

£ = -ei"'\d,.a.x + a^j^ (45) 

TT 

= aoif + -e^'dibj) + -boe'^dittj + -e'^cubj - aif (46) 

TT TT TT 
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There are no time derivatives acting on ao and 60, so they can be regarded as 
Lagrange multipHers. As described in the introduction, they enforce constraints 
on the spatial gauge fields. Although in this case, the constraints look more 
complicated: 

e'-^dibj = -njo 

e^^d.aj = (47) 
The second constraint is easily solved using Hodge's theorem 

a, = d,A+^ (48) 

For hi there will be an additional term in the solution. 

h = 9^n + ^^^ne^^^f (49) 

where A and are periodic functions on the torus and a^, as well as bi, is spa- 
tially constant. So far the analysis followed the standard BF-theory, described 
in the introduction. The new feature is the inhomogeneous part in the differ- 
ential equation for bi. In general, you can always solve these equations with 
help of Green's function, in this case the Green's function for a torus, denoted 
by: The interesting (topological) behavior is given by and 6j. As e"''= is 
invariant under the transformation 

r — > r + L 
L=(L^,0) and £==(0,Ly) (50) 

ips is only allowed to change up to an additional 2tt. With the definition of 

Hi — ^difs the Wilson loops Ai = e'^ ~ and Bi = ~ can be 1 or —1. If 
a quasiparticle encircles a flux quantum, it picks up an additional minus sign. 
Otherwise, it's wave function remains unchanged under a transformation along 
a closed path. 

After inserting this into (I45p . the Lagrangian is simplified considerably: 

C = e^^~^^-Ai7rd,f-d,f)-~^T (51) 

Iji Lij Lii 

Also A can be regarded as a Lagrange multiplier which ensures the current 
conservation, and ^ form a canonical conjugated pair 

[7^,/] = * (52) 

In contrast to the s-wave case, the BF-Hamiltonian is not identically zero. In- 
stead, the conserved current, in this case the spin current, couples to a (spatially) 
constant gauge field. 

n ^ (53) 
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So the Wilson loops, set up by the j^, commute with the Hamiltonian for both 

i, while the ones for the |^ do not. Ai — exp(i <f j^) measures the magnetic flux. 
In other words, the number operator of the magnetic flux can be simultaneously 
diagonalized with the energy. Thus, the ground state is labeled by the number 
of fluxes. Creating quasiparticles does not require a finite energy, therefore 
their particle number is not fixed in the ground state. The total electric flux is 
measured by Bi — exp(i <f j^. Due to the commutation relations of the with 
the Hamiltonian, the ground states can not be labelled by the number of electric 
fluxes. Note, however, that the Wilson loops cannot distinguish a positive from 
a negative flux quantum. At takes the same value for at = n and ai = — tt. In 
addition, you can't count flux quanta with help of the Wilson loops. Two flux 
quanta give rise to a phase of 27r and are, therefore, not distinguishable by any 
other, even number of vortices. Each Wilson loop has the eigenvalues I and -1, 
for even respectively odd number of flux quanta present. 

3.3 Dirac Hamiltonian 

I use (|53p to write down the final result for the Hamiltonian. Note that Oi is 
spatially constant. Thus, the Hamiltonian can be treated as in the field-free 
case. The topological phase only shifts the momentum by a constant amount. 
The Dirac Hamiltonian is given by: 

It is block diagonal and, therefore, splits into two two-component equations. 
The two equations differ only by an overall minus sign, so I will concentrate 
on one. The cigenfunctions of the Hamiltonian can be expanded in momentum 
space by writing 

where are solutions of the Dirac equation for positive resp. for negative 
energies. Here E will always denote a positive quantity. There are nontrivial 
solutions for iJ-l^E)'^ j = if the determinant of the matrix equation vanishes: 

=>E^ ^ vUkx^~^f+vl{ky^^f (56) 

L/'x L/y 

for node I and 3, I denote this pair by (1/3). For the other pair of nodes, (2/4), 
you obtain the same expression but with x and y interchanged. So 




= (55) 
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^(1/3) = JvUk:r-^y+vi(h-^r 



E(2/,) = ^^Uky-fr+vl{k.-f-r (57) 

The total energy is given by the sum of the partial energies. 

As we're on the torus, we need to specify boundary conditions for the spinor 
fields. On the torus, the points x = {xi,X2) and x' = {xi + nL^, X2 + mLy) are 
identified. In general, the spinor fields can difii'er by a phase a at both points. 
The phase of the adjoint fields is then given by —a. Both combined give the 
relation of the spinors at both points: 



*,(f) = U^j{x) (58) 

With U given by: U = diag(e*", e~*", e'", e^*") The Lagrangian has to be the 
same for x and x'. Therefore, the matrices have to fulfill following constraints: 

U^2U = 72 (59) 

The first one is trivially satisfied, while the second one gives a non-trivial con- 
straint on a: 

This determines the boundary conditions to be either periodic or antiperi- 
odic. Which one is chosen is not important as long as I also include all the 
flux sectors. Inserting a magnetic flux changes the boundary conditions from 
periodic to antiperiodic ones and vice versa and will only change the labelling of 
the different sectors but not the overall structure. Here I chose the antiperiodic 
one which sets k to: 

2n + l 

ki = TT 

This choice makes sure that the ground state is unique for ^ = = 0. 

I introduce TTj = fcj — -f^ in order to shorten the notation and to make it 
uniform in the different sectors. The solutions for the spinors u^. and depend 
on the value of k and the Chern-Simons potentials. The cases, where one or both 

of the TTi's are zero, have to be treated separately. However, it is still possible to 
write down a closed expression for the spinor fields with momentum-dependent 
coefficients. 

In the sector — Ay = 1, and also for an excited state in any other sector, 
both TTx and TTy are nonzero. The additional "nodon" symmetries, see Appendix 
B, allow me to treat the two pairs of nodes separately. I will therefore only 
solve the equations for the pair (1/3). Solving for the other pairs is completely 
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analogous. As the theory is not relativistic I chose to normaUze the spinors to 
1 instead of the more common ^2E{k). 

= c{-VATry,E - VFTrx,0,0) 
vl'' = C{E - VFTrx,VATTy, 0,0) 
C(0, 0, - WfTTx, t^ATTy) 

c{0,0,-VA'n-y,E - vpn^) (60) 
with 



2t 

H = 

„2t _ 



^2E{E - VFTTx) 



(61) 



If TTx = and TTj, ^ 0, the energy will be given solely by the canonical 
momentum in the y-direction. In this case, the quasiparticles will be an equal 
mixture of particles and holes. 



= :^{E,7ry,0,0) 

In the case where ttx^O and % = the gap function does not have any effect. 
The Hamiltonian is already diagonal. However, I choose a representation of the 
eigenvectors, which simplifies later computations. 



= ■^{-{E + T,^),E- 1,^,0,0) 

Vl^ = -^{E-TTx,E + TTx, 0,0) 

ul^ = ^{0,0,E-7rx,E + nx) 

vl^ = ^{0,0,-{E + TTx),E-nx) (63) 

Finally consider the zero mode: tIx = = Tiy. In this case, the Hamiltonian 
for the ground state is equal to zero. Thus, the eigenspinors can be chosen 
freely. The simplest choice is most convenient in this case, as it gives the most 
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convenient representation for the spin operators. 



(1,0,0,0) 
(0,1,0,0) 
(0,0,1,0) 

(0,0,0,1) (64) 

In general, the new operators, ai and bi, are a mixture of particles and 
holes, similar to the s-wavc case. In order to express the Hamiltonian in the 
new creation and annihilation, you have to rotate the spinors by a (momentum- 
dependent) unitary matrix U(fc). You have to do this in each sector separately. 
After this rotation, the Hamiltonian is diagonal and can be written as: 

^ = 12^k{('t{k)as{k) + bt{k)bs{k)} (65) 

In all the sectors, apart from the one where Ax = Ay = —1, the energy is 
nonzero for all momenta. Thus, by filling up all the states with negative energy 
and leaving all the other ones empty, one obtains a unique ground state. In 
the forth case, the ground state has zero energy and at first it seems arbitrary, 
whether to fill it or not. For each node, the zero mode can be empty, completely 
filled or occupied by either a spin-up or a spin-down particle. For the node pair 
(1-3) this means a 16-fold degeneracy of the ground state. In order to label the 
ground states one needs to find a set of commuting operators. The eigenvalues 
of those operators don't change under adiabatically switching on a magnetic 
fiux. Later, we will see that only two of the 16 ground states can be connected 
to the ground states of the other sectors by adiabatic processes. 

In a Cooper pair, the electrons have opposite spin and momentum. Breaking 
up a pair and, thus, creating a pair of quasiparticles does not conserve the charge 
but it does conserve spin and momentum. However, the spinors are defined in 
a way that spin-up particles have positive momenta, while spin-down particles 
have negative ones. So momentum does not provide additional information for 
the states and I will ignore it. An additional quantum number is needed to give 
a complete characterization of the states. As mentioned earlier, the chirality is 
an appropriate choice. In the following section I will derive how to obtain the 
quantum numbers and discuss their physical meaning. 

3.4 Symmetries 

The original Hamiltonian is invariant under SU(2) spin rotations, so the spin 
current is conserved and the spin in the z-direction, together with the total 
spin, are good quantum numbers. However, when I introduced the spinors, I 
obscured this symmetry. For a detailed analysis of how to obtain the spin in 



/II ' 



vJ = 
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the different directions, see Appendix A. 



]{al{q)aM-bliq)bsm (66) 
The momentum is also diagonal in the new set of operators. 

= J2ial{q)aM-bl{q)bM) (67) 

s,q 

However, spin-up is coupled to positive momenta and spin-down to negative. 
Therefore, the momentum does not provide any new information and I will 
omit it. 

Analogous to the 3+1 dimensional case, one can construct a fourth 7-matrix, 

73, 

^= = ( -1. 1 ) ' 

which fulfills the Clifford algebra. In the 3+1 dimensional case, the chirality is 
defined by a symmetry transformation, which leaves the massless Dirac equation 
invariant. Its (hermitian) operator is given by the product of all 7-matrices: 




75 = 70717273 
/ 

-/ 



(69) 



In 2+1 dimensions, 73 is redundant for the Dirac equation and can, therefore, 
be combined with 75 to give a bigger internal symmetry group. Both matri- 
ces anticommutc with the other 7-matriccs. Therefore, there is a set of three 
hermitian matrices, which anticommute with each other but commute with the 
matrices ai and a2- Explicitly, they are given by ^73, 75 and 735 = 7375. One 
of these can be included into the set of commuting operators. I will use 75 as it 
is already diagonal. Using another notation, 

Ti = ^7375 

1 

r3 = ^75 (70) 



31 



it is easy to see that they fulfill the spin algebra: 

[ri,r2] = iVs (71) 

The eigenvalue of Fa, I call chiral charge, in analogy to spin and electric charge. 

= al(fc)ai(fc) - 4{k)a2{k) - h\{k)hi{k) + h\{k)h2{k) (72) 

k 

To gain further insight into the chirality, I rewrite it in the old variables. 

= '^It'^I'T ~ ■^34'^3,i - V'il,t'03,T + V'l.^V'i.i (73) 

k 

From ((73)) it is apparent that 2F3 does nothing else than count the difference 
of the number of particles at the node pair (1-3). In contrast, the ordinary 
number operator N is not a good quantum number and does not commute with 
the Hamiltonian. This is what you would expect, since the quasiparticles are 
massless around the nodes. Thus, their number is not conserved. However, these 
quasiparticles are created by destroying a Cooper pair. Creating a quasiparticle 
at eg. node 1 is always connected with creating one at node 3. The difference 
of the number of particles at node 1 and node 3 is unchanged by this process. 

In the isotropic case, vp = va, the chirality forms an U(4) symmetry group. 
For a more detailed analysis of this symmetry see [1] . By letting both velocities 
differ from each other, U(4) breaks down to U(2)xU(2). For each pair of nodes, 
the chirality behaves as the well-known spin algebra. Measurements on cuprates 
show, that the Fermi velocity is approximately one order of magnitude bigger 
than the velocity tangential to the Fermi surface, so the U(2)xU(2) description is 
appropriate. Thus, it is legitimate to regard the chirality as an isospin variable. 

3.5 Ground State Degeneracy 

I will now discuss the ground state degeneracy. The only non-trivial case is 
the fourth sector, where both and Uy are non-zero. In the other sectors, the 
energy gap is always finite. The ground state is determined by filling all the 
negative energy states and leaving the positive ones empty. However, in the 
fourth sector there are states with zero energy and there is an arbitrariness, 
whether to fill these states or not. 

Let fcp be the momentum of this state, fulfilling 

\ ay J 

and let |0) be the state, where the zero mode is empty, 

a,(fco)|0) =6,(fco)|0) -0 
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Creating particles with momentum kg does not require energy. Therefore, acting 
on |0) with a|(fco), ^K^o) or a product of these operators, will also generate a 
ground state. As these operators fulfill anti-commutation relations, you can 
construct 16 linearly independent ground states and label these states with spin 
and chirality. Both symmetries form one triplet, four doublets and five singlets. 
For a detailed analysis, see Appendix C. 

The ground states of the different sectors can be connected by tunneling 
processes as described in [T| . Consider creating a vortex-antivortex pair. Then 
move the vortex along a non-contractible, closed path and let them recombine 
afterward. Moving the vortex around a hole in the torus, creates a magnetic flux 
in that hole and changes the singular part of the corresponding gauge field, 
or (iy, by tt. Thus, doing this adiabatically, connects the different sectors. This 
adiabatic process corresponds to slowly switching on the flux, respectively 
ay. 

For each value of the flux, the spin in z-dircction Sz, the total spin S^, 
the chirality T3, which describes the difference of numbers of particles and 
commute with the Hamiltonian. During an adiabatic process, the states remain 
eigenstates of these operators. 

So first I want to compute the eigenvalues of spin and chirality for the unique 
vacua in the other sectors. As a^lO) = 6^10) = 0, it is straightforward to see that 
the unique vacua are eigenstates to all the operators with eigenvalue zero. In the 
degenerate sector, there are only two states which also fulfill these requirements: 
the first is a linear combination of the completely empty and the completely filled 
state, the other one is a combination of 2 excitations each: 

^1 = (^l + -^aI(A:o)4(M&I(fco)4(fco)) |0) 



1 



(at(fco)6l(fco) + 4(fco)&l(^o)) |0) (74) 



Thus, if you consider the system under adiabatic processes, the huge degen- 
eracy is reduced from 16 to two. Depending on the order in which the fluxes 
were turned on, one of those two is obtained. If the flux in the x-direction is 
switched on first, one obtains ipi at the nodal pair (1,3) and 02 at nodal pair 
(2,4). Taking both nodes into account, there are four states which are singlets in 
spin and chirality. However, only two of those can be reached by the tunneling 
processes I considered. 



3.6 Summary of the first part 

In the first part of my thesis, I derived the effective low-energy theory of d- 
wave superconductors. As it could be expected, the appropriate description 
is a massless Dirac Lagrangian, given by equation (j43p . Furthermore, we saw 
that the conserved current is nothing else than the spin current which couples 
to a topological gauge field = ~^d^j,ips- We then proceeded by analyzing 
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the topological theory on the toru^ and solving for the energy. The results we 
found, differ mainly in two aspects from the ones you expect in a conventional 
topological theory (with gapped quasiparticles) : the ground states in the four 
different sectors are not degenerate, in contrast to the gapped case. In addition, 
in one sector the ground state is not unique but there is instead a rather huge 
degeneracy, in total given by 256 states. Finally, we argued why you could 
focus on only two of these ground states. These were the ones, which could be 
connected to the ground states of the other sectors by tunneling processes. The 
structure we obtained for d-wave superconductor differs from the one for s-wave 
superconductors. It does not show topological features, such as the ground state 
degeneracy on the torus. The second part of this thesis deals with the question 
whether there are other regularities in the theory due to the topology. 



^including magnetic flux tubes through the holes 
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4 Effects of a local perturbation 



4.1 general remarks 

In the first part, I derived a topological description for the low-energy exci- 
tations and solved for the energy spectrum on a torus in the different flux 
sectors. I proceeded by discussing the novel features that arise because of the 
gapless excitations, emphasizing that there are several important differences to 
the fully gapped case. The most significant being, that gapless modes destroy 
the (topological) ground state degeneracy and make the system sensitive to local 
perturbations. However, I will argue that perhaps there are topological features 
left, even in the presence of gapless particles. 

For s-wave superconductors, one obtained a four- fold ground state degener- 
acy on the torus, which was solely determined by the topology [T] . The ground 
states of the different sectors can be connected by adiabatic tunneling processes. 
The energy gap of the ground states in the different sectors is proportional to the 
exponential factor e~^'^"^, where e is the energy of the quasiparticle, n the wind- 
ing number, L the size of the torus and c a normalization constant . Hence, the 
ground states become degenerate in the limit of a large torufl This argument 
does not hold if there are gapless excitations. The quasiparticles are massless 
and free to tunnel, even if the torus becomes very large. Therefore, it is to be 
expected that the degeneracy between the different sectors is lifted. The gapless 
modes also cause one of the sectors to behave exceptionally as they give rise to 
a zero mode and, thus, to a large ground state degeneracy. I argued that only 
two of the states can be connected to the other sectors via adiabatic processes. 
In the following discussion I focus on these two states and their behavior in 
relation to the unique ground states in the other sectors. 

I have already argued that some characteristics of topological orders (such as 
the topological degeneracy on the torus) are destroyed by the gapless excitations. 
However, the superconductor is still described by a topological theory and even 
in the presence of gapless modes, there may be topological features left. In order 
to analyze this possibility, I consider local perturbations and study their effects 
on the ground state energies. In particular, we look for quantities which show 
some regularities in the different sectors. The idea is, to include a magnetic 
point defect at point xq with strength Vq: V = Vq S{x — xq). The coordinates 
of the impurity are arbitrary. They do not affect physical results as the the 
torus is translational invariant. In the fully gapped case, such an impurity 
does not effect the ground state degeneracy. In order to distinguish between 
the different sectors, the quasiparticles have to travel around large cycles and 
measure the flux through the holes. Those paths are not suppressed any longer in 
the gapless case and the local perturbation will shift the energy levels. However, 
it is not obvious, whether these shifts are distinct for each sector or whether 
some regularity reflects the topological degeneracy of the gapped system. 

^In this limit, tunneling processes are suppressed as the action of the quasiparticle S = 
— nL gives an exponential fall-off in Euclidean space 
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The perturbed Hamiltonian is given by: 



n = Hq + V 



J (fx'^i{x) {-iD^Vf^i - iDyV 

+§2(x) {-iDyVfJi - iDxVA^2) *2(x) 




(75) 



The characteristic quantity for the vacuum is given by the Casimir energy. 
It is distinct for each flux sector and is, in general, shifted by the perturbation. 
The shift, however, could be a quantity which is independent of the sectors and 
only determined by the topology. In the following sections I calculate the energy 
shift in the Casimir energy in all four sectors and compare the results. 

The most convenient method to calculate the Casimir energy is, as mentioned 
earlier, by using Green's functions. In section 4.2, I derive how the Green's 
function is related to the vacuum energy: 



The reader may note that special care is needed, when zero modes are con- 
cerned. Whether this mode is filled or not, affects the shift of all the energy 
levels. This contribution can be computed in two ways. One possible way is 
to include this mode in the Green's function. An occupied zero mode simply 
shifts the Green's function by a constant. Another option is to exclude the zero 
mode explicitly and calculate its contribution separately. The Green's function 
representation I use to compute the Casimir energy, is not suitable when there 
is a zero mode present. To first order, both approaches are equivalent and I will 
compute the shift in both ways. However, to compute the shift to second order 
you have to exclude the zero mode from the Green's function as it gives raise 
to a (fake) infrared divergence. 

At this point some general remarks on the interaction potential are in order. 
Choosing a point defect is certainly the simplest perturbation, one can consider. 
However, there are some well-known difficulties with a singular interaction. It 
scatters into all momentum states with the same strength. Therefore, it does 
not provide a cutoff for high momenta and leads to ultraviolet divergencies. 
Normally, this divergence occurs already in the integrals for the first order shift. 
In our case, however, only scattering at the impurity can change momentum, 
whereas the insertion of iujjo at x is momentum conserving. Therefore, to first 
order only scattering processes contribute, where the interaction does not change 
the momentum eigenvalue of the quasiparticles. 

The Casimir energy is obtained by regularizing the vacuum energy on the 
torus as well as in free space. When the vacuum energy on the torus is expressed 
as a sum over windings, the free space vacuum energy is subtracted by omitting 
the term with no windings. As soon as windings are included, the size of the 



E 




(76) 



where G{x, x, w) is the full propagator. 
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torus yields a natural length scale and provides a cutoff for high momenta. This 
is sufficient to ensure a finite expression for the energy up to the first order 
shift, but fails for second order. When the particle scatters two times at the 
impurity, it may have any momentum between the two scattering events. This 
scattering process yields an infinite result, even after the subtraction of the free 
space vacuum energy. The energy shift can be regularized by splitting the point 
of interaction and assuming a finite distance between both scattering events. If 
the quasiparticle does not wind around one of the holes of the torus, then one 
encounters a short distance divergence in the limit where both points coincide. 
However, windings between the scattering events are not sufficient to render the 
integrals finite. The integration over the torus surface diverges also when the 
interaction points coincide. In addition to subtracting the vacuum energy of the 
free space, one must also renormalize the interaction strength by adding local 
counterterms. 

As mentioned earlier, the delta potential does not conserve momentum. 

More specifically, it introduces an interaction between particles at arbitrary 
momentum. Thus, the decoupling of the node pairs is not valid any more and 
particles get scattered between different nodal pairs. To first order, this process 
is irrelevant as iuj^o is momentum independent and the trace ensures an overall 
momentum conservation. There the different nodal pairs can still be regarded 
as independent and the total shift is given by the sum of two independent con- 
tributions. In the second order approximation, the nodal pairs (1,3) and (2,4) 
can mix. As there are two scattering events, the overall momentum conservation 
does not prohibit the particle to be scattered to another node pair. The way 
to implement such processes is straightforward. The divergent behaviour and 
the corresponding regularization techniques are the same for both nodal pairs. 
They differ only in the coupling of the x- and y-direction to the characteristic 
velocities and 7 matrices: vpji and f a72- Therefore, the only difference is given 
in the dependence on the parameters Lx,Ly, vp and v^. 

4.2 Casimir energy 

In this section, I first give a general derivation of how the Green's function (in 
Euclidean space) can be related to the vacuum energy density. We obtain the 
Casimir energy density by subtracting the free space energy density from the 
one on the torus. I also give a derivation how the Green's function looks like in 
our specific case of 2d massless fermions on the torus. For each sector without 
a zero mode, there will be a unique Green's function which will, in general, 
lead to different energy densities in the different sectors. I also derive, how the 
zero mode affects the Green's function and specify it for the degenerate ground 
states in that sector. 

In general, the fermionic Green's function is written as the time-ordered 
product of the spinor fields: 

Ga,;3(x,y,T-T') = (0|T{*„(x,r)*^(y,T')}|0) (77) 

expressed in imaginary time r = it. It is connected to the real-time Green's 
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function by Ga,/3{x,y,T)\^^.^ — iGa,p{x,y,t). The spinors can be expanded in 
the (positive and negative) energy eigensolutions of the total Hamiltonian: 

*„(x,r) = ^*+Jf)e(-4-)+vl/~„(f)e(-.--) 

n 

= -^^a.(p>,(p)e-'?^-~ + 6l(p)i;.(p)e^P~ (78) 
with the convention that 

a,(p)|0) =fe,(p)|0) =0 

(0|at(p) = (0|6t(p)=0 

To calculate the vacuum energy, all the negative energies e~ have to be 
summed up. For that purpose, I assum^ll t < t' and take the derivative of the 
Green's function with respect to t: 

drG^,p{x,y,r~T') = E^"(0|i*Xn(y)*^,™(^)}|0)^^"'"'^'"''^79) 

t<t' — 
n.m 

By using the anticommutation relations it is straightforward to see that 

dr f dxtr{joG{x,x,T-r'))\^,^^ - ^ e^"^ (80) 

n 

This sunjl is still regulated by r' — t and finite as long as r' — r > 0. In the 
limit T t' , the vacuum energy is obtained: 



(fxTi (dTG{x,x,T)jo'j 

d^x f ^iLuTr(^G{x,x,ujh„) (81) 



E 



This integral diverges in the limit t — > 0^ and it must be regularized. In 
addition, one has to subtract the part which corresponds to free space. For 
that purpose, it is helpful to rewrite the propagator. A formal derivation of the 
Poisson resummation formula can be found in Appendix D. However, you can 
also give an argument on how to rewrite the propagator by looking at the path 
integral formalism. 

Consider the probability amplitude K{xa, Xb) for a particle to go from some 
initial point Xa to a final point Xb- 

K{xa,Xb) = e'^'''' (82) 

paths p 



^to choose r — t' < ensures, that the sum runs over the filled levels in the Fermi sea, for 
r > r one sums instead over all empty states 

*In the formula above, you can already see that the Green's function formalism may be 
problematic in presence of the zero modes as it is ambigous whether you include the zero 
mode (and its shift) in the sum or not. 
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As you can only measure absolute values of the probability amplitude, you 
can also add a constant phase to the action S[p]. On the torus, all paths 
can be divided into classes. Two paths are in the same class, if they can be 
continuously transformed into each other. As paths around the handles of the 
torus can not be contracted to a point, it is clear that the classes can be labeled 
by the winding numbers. Thus, the sum over all paths can be broken into two 
sums, where the one over windings encodes all the information of the topology. 
The remaining sum yields nothing else but the probability amplitude in free 
space. In each of the classes, the phase has to be constant, but as the paths 
of different classes cannot be transformed continuously into each other, they 
may be different constants in different classes. In this example, the phase can 
be related to the magnetic flux quanta enclosed by the path hy (f> — i J dPx a. 

In my case, there is also a further contribution originating from the boundary 
conditions. In order to obtain a unique ground state in the non-flux sector, I 
chose the boundary conditions on the spinor fields to be antiperiodic. This 
boundary conditions has the same effect as inserting one flux quantum into 
each one of the holes. More specifically, by inserting flux quanta, one can switch 
between periodic and antiperiodic boundary conditions. Windings around a hole 
with no flux quanta (more generally, an even number of flux quanta) have an 
alternating factor. Adding one flux quanta, consequently, eliminates that factor. 

Now we can insert the (total) Green's function into ([5T|) and compute the 
Casimir energy and its shift in presence of a perturbing potential. The interac- 
tion will be regarded as small and we switch to the interaction picture: 

Hiir) = J d^xe"°^V{x)e-"°^ (83) 

Then, the Green's function in expression (|80p is expanded with regard to 
the interaction strength Vq: 

G{x,x,t,t') = ^(0|T{vI/,,„(f,r)'l',,^(x,T')e-/''"''^^("°^}|0) 

1=1,2 

= ^(O|r{vI/,,,(f,r)^',,0(f,T')}|O) 

1=1,2 

- f dTo{0\T{^,,a.ix,T)^,,pix,T')HjiTo)}\0) 

1=1,2 

J2 j dTl j dT2(0|T{^',,„(x,T)*,,;3(f,r')-ff/(Tl)iJ/(T2)}|0) 



1=1.2 • 



+ 0{V^') (84) 

By using Wick's theorem, one can write the time ordered product as the nor- 
mal ordered product and the sum over all contractions. The contractions are 
nothing but the non- interacting propagators. In the sum, I neglect all contri- 
butions which correspond to disconnected Feynman diagrams. Note that the 
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non- interacting propagator is always "located" at one node paifl I denote the 
pair (1,3) with a superscript (I) and (2,4) with a superscript (II) respectively. 

The impurity is taken to be rigid. It can transfer momentum but not energy 
to the scattering particles. All Green's functions are, therefore, evaluated at 
the same frequency w. This can also be shown formally, by doing a Fourier 
transform of each Green's function. The time integrals then enforce energy 
conservation. 



i=i,ii 

- FoG«(a',a'o,^)7oG«(a'o,a;,c^) 

+ Vi (2^' , a;)7oG(^) (xq, fo, c^)7oG(') (f .S,uj)} 



+ 0{V^') (85) 

With G^^^x, xo,u!) I denote the unperturbed propagator on the torus, whereas 
G{x, y, uj) is the full interacting propagator. For simplicity, I chose to omit the 
index for the different sectors. However, I keep the index for the nodal pairs 
explicit as it will be needed in the calculations for the second order shift. 

To keep track of the various terms in the expansion, it is convenient to 
introduce Feynman diagrams to label the different summands. The notation 
originates from the standard notation for the expansion of a propagator going 
from an initial point if to a final point y. All the propagators are assumed to 
be at the same frequency w. 




Figure 3: standard representation for an interacting greens function, with dis- 
tinct initial and final point 



According to (j8ip . the initial and final point arc identified with each other. 
Thus, each diagram is closed into a loop. Moreover, the operator iuj^q is inserted 
at point X. The trace and both integrations are implicit in all expressions. For 
the Feynman diagrams of the total Casimir energy, see figure S] 

Each propagator G*-'-* (x, y, w) can be written as a sum over windings of the 
free propagators Gq\x + fh^y^uj). The free propagator itself can be separated 
into two pieces, which differ in their short distance behaviour. One couples to 
7o, the other to (a? — y)^. As the commutation relations for the Dirac matrices 
give constraints on the diagrams, it is convenient to label each propagator by 
the Dirac matrix which is picked up. In addition, the winding number can be 
specified. The index for the node pair is, however, kept implicit in all diagrams. 

^The nodal pairs (1,3) and (2,4) are independent without the perturbing potential 
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Figure 4: diagrammatic representation of the vacuum, or Casimir energy as an 
interacting Green's function, where initial and final point coincide 

The free propagator at the different nodes differ only in its dependence on the 
parameters L^, Ly, vp and va- 



- if) rAi) 



i-K^;\\u:\\x'-m (86) 



where x\ = ^ and 2^2 = foi' the nodal pair (1,3). For the other nodal pair, 
the velocities have to be switched. The full propagator is written as a sum over 
windings of the free propagators. 

G«(a-,y,c.) = X^/,;(n)G«(f+n,y,c.) 

n 

= X! 2t:vfvi:^ {-iujK^^'' {\jj\\x + n - y\) 

ft 

\x ~t- y \ 

where the winding numbers are given in multiples of In this thesis, I will 
mostly write x instead of x' and keep the dependence on the nodes implicit in 
the expressions. 

Special care is needed in the zero mode sector. In general, the time ordered 
product of operators is given by the normal ordered product and the contrac- 
tions. When there is a zero mode, this mode can be partly filled. This will 
of course change the notion of normal ordering. When one writes all creation 
operators to the right of the destruction operators, this order will, in general, 
not destroy the ground state if this state is not empty. The zero mode is special 
in that respect, that it does not depend on time or the spatial variables. Thus, 
the Green's functions for the different ground states differ only by a constant. 

G%{x,y,T) = (0|T{vI.,„(f,r)^',^(y,0)}|0) 

= G«(f,2;,r) 

+ (0|vl/- (x,r)^'+ (y,0) - ^-^(y)«'+ (x,t)|6) 

= G«(f,y,r) + Mg (88) 
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The constant contribution becomes simple, when we look at G70 instead as 
M(')7o is diagonal. Using the anticummutation relations, it is easy to see, which 
matrices you pick up when the zero mode is occupied by a quasiparticle: 





E22 


b2{ko) - 


E44 


ai{ko) - 


~E- 


02(^0) - 


-E- 



(89) 

Analogously, for states with several quasiparticles the total contribution is ob- 
tained by the sum of the single ones. For example, the ground states 0i and (1)2 
give the same contribution: 

W '^s 



°Vf" M) = ('? " I m 



The different notion of normal ordering does not give any contribution for 
the non-interacting case, as the zero mode contribution is cancelled by the time 
derivative and does not give any contribution to the Casimir energy. However, 
we will have to take it into account in other orders of perturbation theory. 



4.3 Numerical results in zeroth order 

The Casimir energy in absence of a perturbing potential can be written in a 
closed form, given by: 



i,n 

(91) 

where the winding numbers rii are given in multiples of |^ and fj{n) denotes 
the alternating factors which depend on the sector j . 

After inserting the definition of the Green's function and using the trace 
properties of the 7-matrices, it is straightforward to see that only the term 
coupled to 70 survives. Both integrals can then be performed. The one over 
space is trivial and gives just the surface of the torus. The sum is regularized by 
omitting the term without windings and a finite energy is obtained in all four 
sectors. 



VfVA 
L1L2 

VaL2 

vpLi 



VfL2 



VFL2.-3/2 
vaLi 



)-3/2) 



(92) 
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Figure 5: lowest order graph corresponding to figure ([4]) 

First, one can try to draw some conclusions from the analytic expression. 
For example, one can compare the sectors, where there is one flux tube through 
one of the holes. The energy shift of these sectors differs only in the dependence 
on the torus length. You would expect the ground states in these two sectors 
to behave similarly because of symmetries in the original d-wave descripton. In 
fact, if the length scales are the same — Ly, their energy shift is exactly the 
same. However, the more the two length scales differ, the more differs also the 
Casimir energy in both sectors. This is not a coincidence, but originates in the 
way we quantized the momentum in the first place. After expanding around the 
nodes, we rotated the coordinate system and chose to quantize the theory in 
terms of the new coordinates. Thereby, we assigned the different length scales 
to the node pairs in an arbitrary manner. If we had chosen to quantize in 
the old set of coordinates, the energy dependence on the length scale would 
most probably have been symmetric at the two nodal pairs and, thus, also their 
behaviour for varying this scales. Of course the way to quantize the momentum 
should not make a difference in the final result. However, the description would 
be more involved when quantized in terms of the old coordinate system and 
it is not straightforward to see that the argumentation can be generalized to 
different boundary conditions in a simple way. 

For the numerical calculation, I assume the two lengths of the torus to be 
equal to each other: = Ly — L. Using that choice ensures that the two 
nodal pairs act in the same way. In addition, the velocities vp and ua differ 
approximately by one order of magnitude in cuprates,[]. Therefore, I assume 
^ = 10 in the calculations. After inserting these values, the Casimir energy is 
given by: 

n,m 

( , ^ + , ^ 3 ) (93) 

+ (10m)2 ^(10n)2 + m2 

where and Ay are the Wilson loops introduced in section 3.2. The Casimir 
energy can easily be computed in units of vp/L in the four sectors: 



{Ax, Ay) 


(1,1) 


(1,-1) 


(-1,1) 


(-1,-1) 


E 


-99063.3 


1427.87 


1427.87 


12415.1 
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One finds that the the Casimir energy of the first sectors is lowest, while 
the sector with the degenerate ground states has the highest Casimir energy 
The other two sectors with only one flux quantum have the same energy 
which lies approximately in the middle of the other two. The splitting of the 
ground state energies is algebraic in L, AE ^ 1/L instead of the exponential 
splitting in the s-wave case. 



4.4 Perturbation theory to first order 

The first order shift in the Casimir energy corresponds to the physical process 
were a quasiparticle scatters once at the impurity. It is calculated by inserting 
the first order pertubation of ([55)1 into ([5T|l . Using the trace formulas of the 
Dirac matrices, only two terms survive: 



icoy,, 




icoYo 




duj 
2^ 



-lUJ ^ 



27r 



^ |w| ^{x + n — xo)(xo ~ X ~ rh) 
2-K |a; + n — xollxo — X — to| 



iiTidwIlf + n - fo|)-?^i(kll^o - 



However, both expressions are odd in while the integrals extends over all 
values from — cxd to oo. Thus, every contribution from a positive frequency is 
cancelled by an equally big one for the negative counterpart and the integral 
vanishes. Therefore, the energy shift to first order vanishes in the sectors without 
the zero mode. 

Note, that this argument does not apply for a; = 0, so the contribution of 
the zero mode must be computed in an alternative way. As it is just one mode, 
one can use standard degenerate perturbation theory. The matrix elements of 
the perturbing potential are given by: 

l$\V\^ = ^^at(fc)a,(g>t(fcK(g)-6t(^-)5^(fc)„t(fc)„^(^-) 

k,q,s 

+al{k)bl{q)ul{k)vs{q) + b,{k)aM4{k)us{q) (95) 

To first order perturbation theory, the potential does change neither the energy 
nor the momentum of the scattering quasiparticle. This can be seen from the 

^"with no magnetic fluxes present 
^^as both lengths coincide 



44 



Fcynman diagrams as explained above, but it is also a straightforward conse- 
quence of degenerate perturbation theory. In first order, the perturbation scat- 
ters only to other state in the degenerate manifold. However, when restricted to 
ko, the interaction becomes diagonal in the quasiparticle operators and is given 
by: 

V = Vo^aliko)a,{ko)-bl{ko)bs{ko) (96) 

s 

Generally, the perturbing potential lifts the degeneracy. In this case, the 
degeneracy is only partly lifted. The kets are chosen as the spin and chirality 
eigcnstates which are defined in appendix B. In this representation the impurity 
interaction is diagonal and the energy shift is proportional to the difference in 
the number of particles and antiparticles. 

AE^^ = {i\V\i) 

= VoY,ii\Na,s - Nb,Ai) (97) 

s 

As particles have a "charge" +^ and antiparticles — ^, the energy shift is pro- 
portional to the total charge, i.e. the spin eigenvalue of the state. Consequently, 
the ground states which have S's-eigenvalue are not shifted and behave as the 
ground states in the other three sectors. The two states which can be connected 
to the ground states in the other sectors are spin singletts, and therefore, their 
energy shift (to first order) is the same as for the other sectors. 

However, first order perturbation does not distinguish whether the states arc 
singlets in chirality. So the energy levels are unperturbed also for other states 
than the ones, which can be connected to the other three sectors by adiabatic 
tunneling processes. Thus, it is of interest to consider also the shift in second 
order perturbation theory and see whether this shift distinguishes between the 
chirality eigenstates. 

The same calculation can be done using Green's functions. In first order 
perturbation theory, one needs to consider 

Ai;(l) = -Y, j d'^X j dTodr{lo)pa 

i 

(6|T{M'L(x, r)*,^(y, 0)T/o«'L(xo, ro)vI/,,(xo, ro)}|0) 
= -VoTr^ j cFx j dTo j ^ia;e^'^(^-^°)7o4'^(^+n,xo,a;)7o 

(0|7V[*i,(fo,To)*i/3(y,0)]|0) 
= -VoTrY^ J d^x J ^i(^2n6(uj)"foG'o\x + n,xo,oj)joM'^'^ (98) 

m 

The normal ordered product I abbreviate by M: it depends only on the state 
|0). The sum combined with the integral over the torus gives an integral over 
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all space which can be performed easily: 

j<fxKM\x\) ^ ^ 



Sx'^^Ki{\lo\\x\) = (99) 



J J 27r 

- -FoTr(7oM) (100) 

The matrix M can be read off from ([55]) . If the ground state |0) has the same 
amount of particles as antiparticles, the trace vanishes and the energy shift is 
zero. This method gives exactly the same result as when we calculated the 
energy shift with help of degenerate perturbation theory. 



4.5 Perturbation theory to second order 

The Dirac delta function yields infinite energy shifts to second order in pertur- 
bation theory. As the theory is two-dimensional and its energy relation is linear, 
one expects the divergence to be linear. The special form of the Casimir energy 
takes out the linear divergence leaving a logarithmic one. All the divergencies 
I discussed so far arise in every sector and can be cancelled by local countert- 
erms. In the three sectors without a zero mode, a finite expression is obtained 
after introducing the counterterms. The summation over the windings leads 
to a well-defined expression for the Green's function and the integral over the 
frequency lu can be performed. This is not the case, when there is a zero mode. 
The energy shift involves now the product of three Green's functions: 

E = ^ / d?x j duiicu 

i, 3=1,11 

Tr (7oG(') (f , fo, (xo, xo, uhoG^\xo, x, uj)) (101) 




Even though this does not cause any fundamental difficulties, it makes the 
numerical calculations much more involved. Due to the commutation relations 
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of the Dirac matrices, only certain terms in (jlOip survive after performing the 
trace. All the remaining terms are even functions in w and the integration over 
the frequency gives a finite result. So, second order perturbation theory gives 
the first non-vanishing contribution to the Casimir energy shift in all sectors. 
In fact, there are only four kinds of diagrams which have to be computed as all 
the other terms are cancelled by the trace properties. Only two of them give a 
finite contribution to the Casimir energy after the regularization, see figure ([7]): 



icoYo icoyo 




Figure 7: (1) Figure 8: (2) 

The most obvious divergency arises, when the quasiparticles scatter two 
times without winding around the torus in between. The Green's function 
becomes singular when the two points of interaction coincide. 

Another divergence is given by the spatial integral, which may look surprising 
at first sight. It is instructive to take a closer look on that particular divergence 
in the fourth sector, where the spatial integral (together with one of the sums) 
can be performed easily and gives nothing but a delta function in momentum 
space. 

E = Tr(^ /(fe/ 70*^7oe^^^"e^«"(G(p,c.)7o)'G(<f,c.)) (102) 

The regularization must be done for both momentum integrals. Windings 
in between the scattering events regularize only the loop at xq. 

There are also two asymmetric diagrams which could give a contribution to 
the energy shift. It is easiest to look at the sum of both integrals: 




(103) 



The two diagrams can be seen as the same one, but with different directions 
for the a:— integration. In first order perturbation theory, I argued that such 
diagrams should not contribute to the energy shift. I will show that the contri- 
bution of the sum of those two, indeed, vanishes. 
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All the divergencies I diseussed so far arise in every sector and can be can- 
celled by local counterterms. In the three sectors without a zero mode, a finite 
expression is obtained after introducing the counterterms. The summation over 
the windings leads to a well-defined expression for the Green's function and the 
integral over the frequency w can be performed. This is not the case, when there 
is a zero mode. The zero mode sector is special in that respect that the spatial 
integral can be performed, yielding a delta function in momentum space. 

^ = Y.Y.j4^J^J ^Tr(z.7o(G(fc;-)7o)^Gte-))e-e-^^(104) 

n m 

When you consider only the p = Q part in the propagator, you can replace each 
propagator by -^7o- The expression you obtain is ill-defined and diverges at 
the lower limit of the uj integral: J dw-^. Thus, in the fourth sector you need 
a low-energy cut-off for this integral. Such a cut-off is provided by excluding 
the zero mode. Once it is excluded, the integral over the frequency becomes 
finite. The contribution of the zero mode has to be calculated in an alternative 
way. As was the case for the first order shift, the zero mode is insensitive to 
the symmetries in the theory. Therefore, its divergence is linear, instead of 
logarithmic. 

In the Green's function representation as a sum over windings, this diver- 
gence appears to be an infrared divergence. However, there should not be an 
infrared divergence on the torus as it has a finite volume and, therefore, only a 
finite number of zero modes. The contradiction is solved by realizing that the 
representation of the Green's function as a sum over windings is a good one, 
if the summation converges. This is not the case for uj = and therefore, we 
need to use an alternative way to calculate its contribution to the energy shift. 
When the energy shift of the zero mode is expressed in terms of single particle 
perturbation theory, also an interpretation of the divergence is provided: it is 
then easy to see that the divergence occurs for the scattering into high momenta 
and is not an infrared, but an ultraviolet divergence. 

4.5.1 Proof of finiteness 

Even though the divergent integrals originate from the same short distance 
divergence, I am going to use two different approaches to deal with them. During 
this whole paragraph, I will keep the expressions as general as possible, so that 
the discussion is valid for both nodal pairs and in all sectors. 

First, I consider the divergence which arises from the zero mode. As said 
earlier, in order to provide a low-energy cut-off for the integral over the frequency 
u), I exclude the zero mode from two of the Green's functions. In the remaining 
Green's function I assume that the quasiparticle winds at least once around 
the torus. The case without windings is considered later as it is completely 
analogous to the other three sectors. To see that the divergence is removed, it is 
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sufficient to look at the low-energy behaviour of the regularized Green function: 



G{x,xo,uj) = V'{G'o(f + n,fo,^) - [ d^y Go{y + n,xo,Lu)ll05) 

n 

It has to be finite for lo — otherwise the integral over the frequency is divergent. 

As the sum over windings is responsible for the lu divergence, I neglect 
the first terms and focus on winding numbers where \rh\ ^ \x\. For large m 
the difference whether you integrate over the torus or whether you choose a 
specific point on the torus becomes negligible and the two terms cancel. By 
doing a Taylor expansion around to, I show that the subtraction of the zero 
mode effectively regularizes the Green's function. In fact, it is sufficient to show 
that i^'o(^|"^|) is regularized. The properties for /'iri(cj|n|) follow from it. 

The first non-vanishing contribution from the Taylor expansion is second 
order in x. The first derivative, VKo{Lo\y + n\) as well as the mixed term in the 
second order did2K()(Lo\y + n\) vanish because of symmetry reasons. Therefore, 
the first contribution comes from 

{d! + dl)Ko{u;\y + n\) 

= {-^K,iu;\n\)+c.'Koiu;m^^^f^ (106) 
\n\ LxLy 

Both terms are finite in the limit where a; — > 0. Once the zero mode is excluded 
from one of the Green's functions, the integration over u; becomes finite. 

The missing contribution of the zero-mode is computed with help of degen- 
erate perturbation theory. There it is well-known that the Dirac delta function 
gives rise to a (in our case linear) divergence and regularization techniques have 
been developed, jllj . I excluded the zero mode from the Green's function when 
there were windings between the scattering. The term with no windings must 
be omitted in the calculation. Thus, the energy shift for the scattering in free 
space must be subtracted: 

V ^ ek J (27r)2 efc 

k 

where the sum runs over all states with finite energy e^. To show the con- 
vergence of (jl07p . it is most convenient to expand the integrand in a Taylor 
expansion. As the divergence is linear, there could be a subdivergent term, 
which is logarithmic in q. By powercounting, the logarithmic divergence must 
sit in the first derivative in the expansion. However, this term vanishes because 
of the integration. The first non-vanishing term is proportional to ^ and is 

therefore finite. 

Now I want to consider the term when there are no windings between two 
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scattering events: 

/pOO 

n,rh,l 

K^\\luW{x^x^ + nY)K^\\Lo\^ {x^ -x~lf) (108) 

where the sectors are encoded in /(n, m, l). Setting rfi to zero gives an undefined 
expression. When the second scattering site is displaced by a vector 5x, the 
expression becomes finite for to = 0: 

icoy„ 

n,l 

5x 

K^o'\\u;W{x-xo + nr)K^o\\u;\^{So + Sx-x-r)^) 

(109) 

For small separations, the Besselfunction can be replaced by its asymptotic 
behaviour: 

K^o^\\ujy^iSxy\) - H\uj'\) + ln(2) + *(1) - In |W for Sx 

where ^'(1) is the Euler constant and is labeled by the nodal pair, uj' is 

given by uj' = co . The first three terms give finite contributions for the case 
where there are no windings. The last can be cancelled by a local counterterm. 

Ha. = V„H\x-Xo)Mai^.To)adro'f,^ix,T) (110) 

which gives a non vanishing contribution in first order perturbation theory: 

aV„^TiiJ d^x J §^tiujoG<^'\x,xo,iu)io^oG^'\xo,x,u;)) (111) 

The expression above can be separated into two terms. First consider the term 
where both propagators contribute a 70 to the trace. That term cancels exactly 
with p09p if a is chosen to be: 

a = -;^(ln|Jf|(^) +ln|(5f|(^^)) 

The term, where both propagator contribute with a (x — xq)^ is needed for the 
second diverging diagram: 
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Figure 9: Feynman diagram corresponding to Hp 



icoy„ 




n.m.l 



(x — Xo + n) {x — xq~\- I) 



\x — xq + n\\x — xo + I 

K['\\u;Wix-xo + 7ir)K['\\Lo\^{xo-x-l)^) 

(112) 

As in the previous diagram, this diagram has a divergence when the quasiparti- 
cles do not wind between the scattering events. We spht the point of interaction 
in the same way as above: 




= 4 



d^x 



(x-Xp+n) (x-Xq+I) (j) r-^ , -..2 \ 

\x-xo + n\ \x-xo + l\ 



K 



^Iv {xo + 5x - X - ly) 



(113) 



The term where m = is treated in the same way as above. The Besselfunction 
is replaced by its asymptotic behaviour at the origin and the singular piece is 
cancelled by Ha- 

The other divergence arises for m but fi = I — Q. In order to obtain 
the same divergent behaviour in a first order diagram, this counterterm needs 
to couple to 7V, see figure ([9]): 



Hp = V^5^{x-XQ)-^l„{x,TQ)p-fQ^VQ^i„{x,TQ) 



(114) 
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This counterterm includes also finite contributions to the energy shift: 



4/3 / d^x I duo 



(115) 



1^1 
(27r)3 



K^\\u:W{x~^x^ + nY)- 



K't\\ujW{xo^x + 7nY) 



\xa — x + m| 

K^\\u;W{x~xo + nf)\uj\K^\\u:W{x^-x + mY 



+ k'^\\u;W{x~Xo + n)2)|c^|i^W(|a;|v/(fo ~ f + m)^) 



^ (x — xq + n) — + ; 
|x — xo + ri| |x — xo + ' 



Only the first term is logarithmic divergent when to = 0. When (3 is set to 
/3 = — ^ it cancels the divergence in (|113p in the limit where (Jx ^ 0. The 
spatial integration over the sum of both divergent terms becomes finite: 

d^x{Ko{\uj{x~xa + ^)i) -^i(M^-^o)l) _ ifo(|c^ri|)i^i(|cc;(f - fo)l)' < oo 



\u}{x- xo)| 

The divergence arises for small |x — xol, so by reducing the area of integration 
such that K[){\uj{x — xq + n)\) « Xod'^'^l the leading divergence cancels. Note, 
that it is not important at which nodal pair the Green's function is located. 
Their contribution differs by a finite amount, while the divergent part is exactly 
the same for both node pairs. 

Finally, I consider the sum of the asymmetric diagrams: 



icoYo 



icoyo 




K['\\Lum\){K','>{\u;ix~Xo + n)\) 



jx — Xo + n\ \rn\ 



K^'Wcoix-xo + m)|) - i^^*'(|a;(x - .xo + n)\)K['\\u;{x - xq + n)\} (116) 

For TO 7^ the diagrams vanish because of symmetry reasons. For zero wind- 
ings, the diagrams in equation (jll6|) are not properly defined as one of the 
Besselfunctions diverges, while the expression in the parenthesis becomes zero. 
When the point of interaction is splitted, (|116p becomes finite and one can 
study the limiting behaviour for dx 0. For the nodal pair (1,3), dx is given 
by (Sx/vft^u/va)- For the other node pair, the velocities are switched. In 
the following discussion, I also include a subscript on 5x to keep track of the 
different dependence of the velocities: 



52 



it, I 

(117) 

The limiting behaviour for the Bessclfunctions is known, each diagram for itself 
diverges linearly in \Sx\. The sum of both diagrams will take out the leading 
order in the divergence. But there may still be a divergence in subleading order 
left. Therefore, I analyze how the expression in the parenthesis acts for small 
separations \Sx\. For derivating the Bessel functions I use: 

K'oix) = -K,{x) 

K[{x) = -Ko{x)--Ki{x) (118) 
x 

so that the result can be written in terms of these two functions, only. The total 
contribution from the two asymmetric diagrams can be written as: 



Trj:JrfJd.xi^) 



LO 4 \5xi 



27r \5xj\ 

{x-xo+n)kj{x-xo + l)^i {i),. ,^ ^ ^\\\i^{i)n - , r\\\ 

\x-xo+n\ + 

((x-xo+n)Aj(f-fo + n)A,-) ^^(,)^,^^^ ^ ,^M^rX^)/ 
\x — xo + n\- 



+ Hi — , -.|2 —Kq '{\uj{x -xo + l)\)K^ (|<^(a; -xo+ n)\) 



^ {{x -X0 + n)Aj){{x -xo+ n)Ai) (j) , r^t, K['\\u){x - xp + n)\) 

+ 2 \x-xo+n\^ ^° Mx-xo + OI) |^(^_^^+^)| 

- A. . aXHHx-xo + r) ^''?r^l' !° + "^'^ (119) 

|a;-a;o + n| 

Where A is a unit vector with direction Sx. The first two terms give finite 
contributions to the energy shift, but the others exhibit a subleading divergence 
which is logarithmic in x. The counterterms which we have introduced so far are 
not able to cancel this divergence as it is direction dependent. Point splitting 
introduced a specific direction which was not present in the original graph. This 
can be remedied either by integrating over all directions or by a counterterm 
which is also direction dependent. The former is more intuitive. Moreover, it is 
quite straightforward to see that the divergent part vanishes, when integrated 
over all directions. However, it is much simpler to compute the finite contribu- 
tion with help of counterterms and, therefore, I use this method instead. In the 
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end, both methods give the same result, namely that the asymmetric diagrams 
do not contribute to the energy shift: 



(■p — ^ p — . ^1 p — ^i3-)^i (1^(2; - xo + m)|) 



(x-xo+to)7 (i) ^ ^ i^f^(|tj(x - To + ™)|) - fo + m)Ai 

— — _ I (I^^Ka; - To + m)|) — — — ^ ) p — , ^1 } 

|t — To + to| \ijj[x — xq ^ m)\ |t — tq + m| 

?5"5°^? i^f^(Mx-fo+n)|) '"'(!-? + ?^ i^f)(Mf-fo + m)|) 
|f-fo+n| 1 ^' ^ 1 VI V /I/ 



4 



(2^) 



3 



((t - fp + m)Ai)((f - To + m)ii) (,) ^ -m\s-W/i - , -mn 

- , -.12 -^0 - XQ + '(|w(t - To + m)|) 

(t - To + n)e^ (f - To + m) A, (,) ^ -mxt^W/i - , -m^ 

+ —r^ — — ^T— "1 - To + n)\)K\ {\uj\x - xq + to) ) 

|f-To + n| |f-fo + TO| 1 ^1 ^ " nj 1 w \ j\j 

((f-fo + TO)A^)((f-fo + TO)ei) (i) ^ xf^(|w(f - To + to)|) 

- 2 — — _ _ K^'(\uj[x- xo + n)\) — — — _^ ^ } 

|T-To + TOp |w(t — to + to)| 

+ eiA.if^ ^ T - To + n ) -— — ^ } 120) 

\UJ[X — To + TOjl 

In order to cancel the divergence in (|119p . must be chosen as 



ei = - 



- E 

27r ^ 



\Sxi\ 
^ |(Sf,| 



A 



In that case, the countcrtcrm cancels not only the divergent term but also all 
the remaining, finite terms from (jll7[) . A finite contribution fi'om the asym- 
metric diagrams would have been inconsistent with the underlying symmetry. 
Furthermore, there is no specific direction in the model. A local impurity on 
the torus cannot introduce a direction dependent energy shift. 



4.5.2 second order results 

In this section, I am going to summarize the results for the energy shift to second 
order perturbation theory. The analytic expressions I obtained are finite and 
can be computed numerically. As there are a lot of terms to take into account, 
I organize the results by the diagrams which they originate from. The total 
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energy shift to second order in Vq is most conveniently written as a sum of five 
diagrams: 




(121) 



The first two contributions are analogous in all sectors. In the previous section 
I derived that both AEa and AE/j are finite in all sectors. In the three sectors 
without a zero mode, AEr does not need regularization. In the zero mode sector, 
I excluded the zero mode from two of the Green's functions. In that sector, I 
need to compute the shift AE!^, which is due to the zero mode, separately. I 
will use the notion of the Wilson loops and Ay, introduced in section 3.2, to 
give an explicit dependence of the energy shifts on the sectors. 



(HIuWV/vfva) - ln(2) - *(1)) • 
{K^o'\\co{x~xo + n)\)Kl^\\uj{xo - x - l)\) 

{x-Xo + n){x-Xo + l) rj) _ -MM^(j)/| - TmM 

|a;-a;o + n| + 

(122) 

AEp = 4jd'xl<hji^r^^'i-A,r^{-Ayr^[ 

J2k!^\\^^\)4'\Hx-xo)\)4'\H^o - x)\) 
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+ {^if«(|c^m|)ifp)(|c.(f-fo)|) 

3 

k['\\u{xo-x)\) 

i n 

^1 (It^ia^o -x + m)|)— , ^1 

- i^^*)(|a.(f-fo+")|)i^,^'^(k(a'o-:?+m)|)}] (123) 

In the non-zero mode sectors, the contribution from the other summands is 
trivial to write down: 

Kl,'\\um\)Klf\\Loix - xo + n)\)K^o'\H^o -x-T)\) 

k['\\uj{x-xo + n)\)K'f\\uj{xo - x - t)]) (124) 

In the zero mode sector, I replace the non-interacting Green's function with the 
one, where the zero mode is excluded: 

G{x,xo,uj) = V'{G'(f + n,fo,^) - ^ I d'^yG{y + n,xo,uj)} 

n 

= E^ Jd'yh'^MKoiHx~xo)\)~Koi\u;y\)) 

n 

and insert the modified Green's functions in p2H) 
AEr = ^jd^xj<fyj<L.{^fY.'T.'K'^' 



{K^,=\\u:{x~x, + n)\) Ki^\\u^{y + n)\)} ■ 
{K^,'\\u^{x, - f - 01) - Ki=\\u{y- m 

+ { i":' t % k[=\W{x-xo + n)|) - p-l.K,{\u:{y + n)\)} 
\x — .To + n\ \y ~\- n\ 
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{ ,:° : ' Ki'\H^o - x-r)\) - l-^K[^\\u;{y- r)|)}](126) 

\xo-x-i\ \y-i\ 

Even though I replaced only two of the three Green's function^, one finds 
that the energy shift of all the filled levels is independent of whether the zero 
mode is (partly) filled or not. However, the energy shift of the zero mode itself 
to the Casimir energy is, in general, distinct for the different ground states: 



V ek J (Stt)^ efc 

k 

The expressions do not show a simple dependence on the windings. It is 
not possible to extract some significant features from the analytic results. In 
addition, the computation of the sums is very time-consuming. For low frequen- 
cies, the convergence of the sums is very slow. In addition, the expressions are 
given by modified Bessel functions. While this does not give any fundamental 
problems, it makes the numerical calculations very slow. With a computational 
program like mathematica or matlab, the expressions can not be computed in 
an acceptable time intervall. 



4.6 Conclusions of the second part 

In the second part of this work, we considered a local perturbation, namely 
a magnetic point defect, and studied its effects on the Casimir energy in the 
different sectors. The total Casimir energy can be written in terms of Green's 
functions and be calculated perturbatively. I showed that the Casimir energy (in 
case of no perturbing potential) is distinct for different sector^ and looked for 
regularities in the energy shift induced by the magnetic defect. In the discussion 
I focused on the three unique ground states in the non-zero mode sectors and the 
two (degenerate) ground states in the zero-mode sector which were connected 
to the other ones by adiabatic tunneling processes. 

A major result of this discussion is that there are, indeed, some regularities 
for these ground states as their energy shift to first order perturbation theory 
vanishes. Note that the first order shift is not generically zero, but that there 
are certain ground states energies which are shifted by the perturbation. 

I also derived the analytic expressions for the energy shift in second order. 
As the singular perturbation introduced new UV divergencies, I renormalized 
the interaction and showed that the contributions are finite in all sectors. To 
write the programs which are necessary to calculate the second order shift is 
beyond the scope of this work and will be left for later work. 



^■^Tho Green's function G{xo,xo,u!) is unchanged and could, in principle, give different 
contributions for the different ground states. 

^''as could be expected in presence of gapless particles 
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5 Final conclusions and Outlook 



The 2+1 Dirac theory proves an appropriate model for the d-wave supercon- 
ductor. In general, a two-dimensional representation is sufhcient to describe a 
two-dimensional Dirac theory. The four-dimensional representation has, how- 
ever, several advantages, when it comes to the interpretation of the physical 
properties. In a two-dimensional representation, there are no closed expressions 
for the spin charge or the spin current. Combining the fermion fields of opposite 
nodes into a four dimensional spinor, provides a simple interpretation for the 
density as the spin and the conserved current as the spin current. 

An equal copy of the representation arises for the other node pair. There 
are various versions of how to distribute the fermion operators into the spinors 
and the suitable choice depends strongly on the quantities you want to study 
and which symmetries you want to keep manifest in the model. As I wanted to 
study the coupling to magnetic flux, it was most convenient to take an approach 
which masks the spin rotational symmetry. On the other hand, it allows us to 
identify 5*^^ with the spin density and to write down the spin current in a 
simple form. 

The most delicate part in the derivation of the Dirac theory, is the point 
where you have to eliminate the topological phase in order to obtain spatially 
independent matrices. As I already mentioned earlier, the naive choice to dis- 
tribute the phase equally to spin-up and spin-down leads to multi-valued spinors. 
To distribute it solely to spin- up particles gives a reasonable result, but one has 
to keep in mind that this induces an artificial asymmetry in the behavior of the 
fermion operators. The spin current is only coupled to the topological phase, 
a result which remains true, even if you choose another way of attaching the 
topological phase to spin, see [1] . 

The effective theory you obtain, looks like an ordinary Dirac theory for mass- 
less spin 1/2 fermions. In order to analyze the theory for topological properties, 
it is convenient to consider a torus as the underlying manifold. Further, you ne- 
glect vortices to simplify the analysis. Magnetic fluxes arc now implemented by 
flux tubes through the holes of the torus and can be measured by taking quasi- 
particles along non-contractible paths around these holes. The gauge field, that 
enters the effective theory, can be regarded as spatially constant and classifies 
four different sectors, one for each configuration of fiux quanta. For an s-wave 
superconductor, the structure of each of the sectors is identical. The ground 
state can be labeled by either the magnetic or the electric fiux as both, quasipar- 
ticles and vortices, are gapped. In the d-wave case, only magnetic fluxes can be 
used to describe the different topological sector. You find, that one of the sec- 
tors behaves exceptional. Which one it is, depends on the boundary conditions 
you impose on the fermion operators. The ground state is highly degenerate 
in this sector. It is 16-fold degenerate for each pair of nodes. The complete 
degeneracy is, thus, given by 256, as both pairs can be treated independently. 

In order to describe the degenerate ground state, the observables spin, mo- 
mentum and charge are not sufficient. As stated earlier, there is also another 
quantity, besides spin and momentum, which is conserved. The chiral charge 
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measures the difference in the number of quasiparticles at opposite nodes. It 
provides the missing tool to label each one of the degenerate states uniquely. 
However, by changing boundary conditions, you can change which of the sectors 
behaves exceptional. This suggests, that this (huge) degeneracy in one of the 
sectors might not be fundamental. As in the s-wave case, you can connect the 
different sectors by tunneling processes. Doing this adiabatically, in particular 
when going from a (unique) ground state to the degenerate sector, you pick out 
one particular ground state out of all degenerate ones. One finds, that the ob- 
tained state depends on the sector you started from. In the end, there are only 
two states which can be obtained by adiabatic tunneling processes. They look 
very similar to each other and can be translated into each other by replacing 
a; y. 

As the quasiparticles are masslcss, the ground states in the different sectors 
are not degenerate and generally the different sectors will mix due to tunneling 
processes. It is an interesting question whether some regularity still reflects 
the topological order. In this thesis, we chose to examine the response of the 
Casimir energy to a local impurity. For gapped quasiparticles, local impurities 
do not affect the ground state energy. We found that this is also true in case of 
massless particles to first order perturbation theory. The Casimir energy of the 
two states which could be connected to the ground states in the other sectors, 
are not shifted and, thus, behave as the ground states of the other sectors. 
In second order perturbation theory. I derived the analytic expression for the 
energy shift in second order. The calculations become much more involved and 
the expressions become divergent. The divergencies can be canceled by local 
counterterms and the interaction strength Vb is renormalized. The expressions 
show no simple dependence on the winding modes and there is no possibility to 
make statements on the second order energy shift without numerical calculation. 
However, the numerical results require a lot of time and computational power 
and will not be presented in this thesis. 

Some experiments indicate that cuprates do not have a pure d-wave pairing 
but rather a mixture of d- and s-wave pairing. As the cuprates are anisotropic, 
the relative strength of the order parameters depends on the direction and is, 
therefore, tunable. In appendix D I show, how the excitation spectrum changes, 
if a (small) s-wave contribution is added. It would also be interesting to see, 
what happens to the degenerate ground states if you (slowly) switch on this 
s-wave contribution. With an s-wave contribution, there is a unique ground 
state in all four sectors. By a perturbative calculation, one should be able to 
find out, which one is selected. 

Another interesting feature is to discuss the effect of shifting the nodes away 
from the position k = (ztkp,ztkp). This certainly effects the energy spectrum. 
The simplest possibility is just to rotate the nodes with a small angle. This 
should not change the structure of the ground state in a qualitative way. Actu- 
ally, by expanding the field operators around the nodes and changing to relative 
coordinates, namely changing from {qx,<ly) to (^'1,(72)5 the absolute position of 
the nodes becomes less important and is reflected only in the parameter i;a- The 
only needed feature is the symmetric arrangement of the nodes. This should 
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lead to a certain robustness of the structure against rotations. Deviations of the 
symmetric arrangement, however, should affect the ground state structure. 

In addition, I completely neglected vortices on the torus by setting the vortex 
current = 0. By allowing vortices only as thin fluxes through the holes, the 
theory remains regular on the surface of the torus. The next step would be to 
consider stationary thin flux tubes on the torus. If the core of the flux tube 
is much smaller than the magnetic length, it is still a good approximation to 
consider the absolute value of the order parameter to be constant. 
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A Spin Current 



As mentioned in the main text, the spin current is the conserved Noether current 
of our model However, our choice of spinors masks the invariance under the 
SU(2) spin rotation group. In order to derive the Noether current, another 
choice of the spinor representation 3J is more useful as it leaves the rotational 
invariance manifest: 

= (4T'4i'^J+2,T,-V'j+2,i) (A-1) 

The minus sign for ^/'j+2.x was chosen so that transforms as a spinor. With 
this convention, spin transformation can be written as matrices, but the chiral 
symmetry is obscured. 

With this choice of the spinor, it is now easy to formulate the operator for 
spin rotations. It is obvious that the Lagrangian will be invariant under this 
(global) transformation. 



U = exp ( ^6* • E,; ) (A-2) 

where 



2 V 

and fji being the well-known Pauli matrices: 



= (? o) 

= ( ° V ) 

- (J -l) (^-3) 

With our choice of spinors, see the spin transformations are more cum- 
bersome to write down because of two reasons: first, we chose the spinors to 
have definite spin, ^E" having spin ^ and ^ having — ^ . Therefore, a spin ro- 
tation mixes ^ and ^E" and thus, it is not possible to represent it as a matrix 
as in [3]. Secondly, although "^{q) looks local in q-space it is not. Note the 
convention, I used to define the fermionic operators ()20p . q was defined in such 
a way, that it gives the direction of the spin current. Spin rotations on the other 
hand are local operations. For example flipping the spin of a particle reverses 
its contribution to the spin current, changing q to —q. 

In order to derive the correct form of the spin operators, it is convenient to 
write down the transformations for the original fermionic operators Co-(fc) and 
translate this expression to the linearized ones. The original operators transform 
in the usual way under rotations with angle a: 

(cT,Ci) ^ exp (j|crj)(ct,cx) (A-4) 
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The spin operators themselves are given by 

k ^ ' 

Using equation (j20p . these expressions can be written in terms of the hn- 
earized fermion operators. After transforming to x-space, one has to attach 
the singular phase to spin- up operators. Then, I use the relations to ([M]) . 
to obtain the final expression for the spin operators. I do this calculation ex- 
plicitely for Sx- Sy is obtained in a completely analogous way, while Sz is even 
simpler to derive. 

Due to the definition of the fermionic field operators, the spin in the x- and y- 
directions looks nonlocal in the linearized version. Although the sign convention 
may seem unnecessary complicated it is nevertheless crucial, to write (jl9p as a 
Dirac equation. 

q 

+ V[x(9)V'i,T(-'?) + ^li('7>3,T(-'7) (A-6) 
and after transforming back to x-space: 



+ ?Aj^(f)V'i,T(-^)+ ^14(^)^3,1 (-^) (A-7) 

The next step is to attach the singular flux to the spin-up field operators. 
I neglect regular gauge transformations, as their phase contributions cancel 
anyway. For clarity, I introduce a new notation for the fermionic field with 
attached flux: 



X,a{x) = e-^'^=(^V,,T(^) (A-8) 
Xjdi^) - V'ja(^) (A-9) 

The Fourier transforms of the ^ fields are defined exactly as for the -0 flelds, 
namely that both spin components transform with opposite momentum. 

x,,t(^) = E^''"xm(9) 



9 
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However, as the singular phase is space dependent, it changes the coupling 
of the fermionic fields. Without a singular gauge transformation, spin-up and 
spin-down couple with "opposite" q values. The presence of the gauge field 
shifts one of these momenta. 



q.k 



+ (y d':re'^(-)e'('~+'=)-^)(xl,x(fc)Xi.T(9l+xLW^^^ (A-11) 
In order to obtain the shift, one needs to find the value of the functions 

f-iip,q,k) = [ (fxe'"^'^'-'^e'"^^+^'>^ (A-12) 



You can take tps {x) to be given by: 

ifs{x) — ~2axX — 2ayy 

After inserting this expression into (|A-12p , the integrals can be computed easily 
to give: 

f+iv, q, k) = /_ (if, q, k) = 5^,_k+2a 

In the following discussion, I keep writing q instead of fc 4- 2a in order to 
shorten the notation. The reader should, however, keep the equivalence in mind. 
Note, that the spin operator is now non-local in momentum space because of 
the attachment of the singular phase. The remaining step is to express the 
obtained expression for the spin in terms of the quasiparticle operators, and 
hi which diagonalize the Hamiltonian. Even though the eigenspinors form an 
orthonormal basis for each momentum, scalar products such as Wj(fc) • Vj{q) are 
generally non-zero. 

First note that Xj,iy{k) and Xj,cr{q) are linear combinations of quasiparticle 
operators with the same energy: 



E{k) = Uikx-^Y+vl{ky~^Y 



ivlhkx + -^r + vi{-ky + -^y 



E{q) (A-13) 
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Assume that neither ky — ^ nor fca; — -^^ is equal to zero. To perform the 
computation, it is favorable to introduce 2-dimcnsional vectors: 



wi{k) = c(-VAiky-Y'),E-VF{k.j;-^) 

W2{k) = c(^E-VF{k,-^lvA{ky-^)^ (A-14) 

with the same normalization constant c as in (|60p. The vectors fulfill the rela- 
tions: 



wi(k)-wi{q) = 
W2{k) ■ W2{q) = 

Wi{k)-W2{q) = - 



1„ "-y 



W,{k)-W,{q) ^ ' (A-15) 

As the spin operator is summed over all momenta, we are allowed to relabel 
k ^ q and it can be written as 

Sx ^ xI,T(^)Xi,i(9)-X3a(^)X3,T(9) + Xu(^)Xi,T(<z)-X3,T(fc)X3,i(9) 



+ (xU(4x3,TW)-(4jt;||) 

E (aI(^K(fc) + fei(^)w2(fc))-(-4(9>i(9l+^i(9>2(g)) 

+ (-a2(^)«;t(fc)+6t(fc)«;|(fc))-(ai(g)u;i(g) + 6^(g>2(g)) (A-16) 
By using the relations (lA-lsp . one obtains a rather simple result: 

= E sign(fcy-^)(4('?)&I(fc)+«I(9)^2(fc)+&I(fc)a2(g) 

ky-^^0 

+ bl{k)aM) (A-f7) 



The same calculations can be done for (p^ and In the case that 

~ TT 7^ 0, you obtain the same result as for ([60| . In the case that only 
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ky — jf- =0 but not kx — jf^, you obtain the same result but with x <-*■ y 
interchanged. All these cases can be summarized by introducing 

/(fc) = sign(fc, - ^) if ky-^^0 

/(fc) = sign(fc, - ^) if ky^^^O (A-18) 

In the last case, I chose the eigenvectors in such a way, that the spin 

operator looks the same for all momenta: 

al{ko)bl{ko) + (fco)?^2(fco) + b{{ko)a2{ko) + 5^(A?o)ai(fco) 

The computation of Sy and Sz can be done in a similar way. The operators 
can be written in a compact way as: 

S. = J2 /(^) {<^UQ}b\{k) + a\{q)b2{k) + b\(k)aM + h\{k)aM) 

+ al{ko)b\{ko) + a\{ko)h2{ko) + b\{ko)a2{ko) + 6^ (A?o)ai(A?o) 
Sy = » E /(^) (-"k^l&Ufc) - a\{q)b2{k) + ?>I(fc)a2(g) + 

k^ko 

- a\{ko)b\{kQ) - a|(fco)fe2(fco) + b\{ko)a2{ko) + 4(^o)ai(A;o) 
5, = ^al(^)ai(fc)+4(fc)a2(fc)-6l(fc)5i(fc)-6^(fc)52(fc) (A-19) 

B Symmetries of the Effective Lagrangian 

The effective Lagrangian exhibits additional symmetries, which were not present 
at the original Hamiltonian. This is mostly due to that the node pairs (1-3) and 
(2-4) are completely decoupled after expanding Cfc and c\. around the nodes. 

Thus the Lagrangian will be invariant if we apply separate SU(2) spin ro- 
tations on the two pairs of nodes ^E"! and ^'2. This is also true for U(l) charge 
transformations. As was pointed out in [3], this implies two new conserved 
charges, called nodon charges, It was also mentioned that additional in- 

teraction terms, e.g. Umklapp scattering processes, will not conserve the nodon 
charge while being consistent with the original U(l) and SU(2) symmetries. 

In our case there is even one more symmetry due to our choice of distributing 
the phase of the gap function A{R). 

V^.a ^ e^^^"/') (B-1) 

(|B-ip introduces a new gauge symmetry, where the singular and the regular 
phase are transformed simultaneously: 

fr Vr - l3 

ips Vs+P (B-2) 
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Table 1: eigenvalues of the degenerate ground states in sector 4 



Written in the spinor fields "ifj this looks like the U(l) gauge transformation. 

e^^/^^-j (B-3) 

Again two different values of P are allowed for the different pairs of nodes. 
The invariance under U(l) charge transformations implies that the ^E" fields are 
neutral while this obviously does not apply to tpj. 

C Analysis of Spin and Chirality 

Here I want to sketch, how to derive the structure in sector four. Due to the zero 
mode, there are 16 linearly independent states with zero energy: these consist 
of quasiparticles with momentum ^0 = (x^, x^) and can be labeled by spin and 
chirality, see Table [TJ 

Both, spin and chirality, can be represented by one triplet, four doublets and 
five singlets. There are four states with eigenvalue zero for spin and chirality. 
In order to see, which ones are singlets and which belong to the triplets, you 
need to compute the eigenvalues to the Casimir operators, and F^: 

= r^+r^ + r^ (c-i) 

Let's look at the spin part first. I rewrite the Casimir operator in terms of 
the ladder operators, S+ and S- 

S+ = {Sx + iSy) 

= E/W(4(9)&2(fc) + 4(g16iW) (C-2) 

k 
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S— — {Sx iSy) 

= 5^/(fc)(6l(fc)a2(^) + 6j(fc)ai(9-)) (C-3) 

k 

Singlets will be eigenstates to with eigenvalue zero. 

S+S-\Q) = 
S+S.a\b\\0) = [a\b\-albl)\0) 

S+S-albl\0) = (albl-a\b\)\0) 
S+S-a\alb\bl\0) = (C-4) 
Thus, three states are singlets in spin, namely: 

|0) 

^(al6l+46^)|0) and 

a\4b\bl\0) (C-5) 

The remaining state, 

i=(4&l-al4)|0) 

forms a spin triplet, together with a\a\ |0) and b\bl\0). 

The same analysis can be done for the chirality operators. They can be 
computed by Fi = ^^''^7375^' and r2 = |4'^73^'. Explicitly, they are given by: 

Ti = ^ ^ (-4(^)02 + 61 62(^)6! (^) + 4 (fc)ai(^)) 

+ I (a\(ko)bl{ko) - 6i(fco)a2(fco) + al{ko)b\iko) - 62(^0)01(^0)) 
r2 = ^ ^ (-4(fc)a2(fc) + 6i(fc)6^(fc) + 62(fc)6l(fc)-4(fc)ai(fc)) 

+ I (a\{ko)bl{kd - b,{k,,)a.2[h,) - al(ko)bl(ko) + 62(^0)01 (fco))C-6) 
The ladder operators are computed analogously to the spin case: 



r+ = -i J2 b\{k)b2{k) + al{k)ai{k) 

k^ka 

+ i (a\{ko)bl{ko) - 61(^0)02(^0)) 
T_ = iJ2 a\{k)a2{k) + bl{k)bi{k) 

k^ko 

+ i (al{ko)b\{ko) - b2{ko)ai{ko)') (C-7) 
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Again, one can act with the ladder operators on the states, computing the 
effect of r_|_r_. They pick another state to be the member of the (chirahty) 
triplet: 

r+r_|o) = {l - alalblbl) \0) 

T+r_a\b\\0) = 
T+T_albl\0) = 
T+r_a\alb\bl\0) = [a\4b\bl -l)\0) (C-8) 
Thus, there are also three chirality singlet states: 

^ (l + a\alb\bl) |0) 
a\b\\0) 

44lo) (c-9) 



The triplet is given by 



4b\ 



^ (l - alalblbl) |0) 

a\bl\0) (C-10) 

Now it is easy to sec the structure of the ground state: there are two states, 
which are singlets, both in chirality and spin: 

^ (l + a\alb\bl) |0) 

^(46t+44)|0) (C-U) 



The state 



(l - a\alb\bl) |0) 



\bl - alb\) |0) 



is singlet in spin but not chirality and the state 

is singlet in chirality but not in spin. These ground states can now be compared 
to the (unique) ground states in the other sectors. 

D The Poisson resummation formula 

The Poisson resummation formula gives a relation between a function / and its 
corresponding Fourier transform /. 

oo oo 

J2 f{x + m)= J2 /Ne^'"'""^ (D-1) 
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where both m and n are integers. This equahty can be used to transfer the sum 
over momenta to a sum over windings on the torus. First, I express the Green's 
function in terms of dimensionless variables : n and x' . 



G{x,y,uj) = Y.Gip,Lo)e^p{ip{x~y)) (D-2) 
p 

G{n, Lu) exp (27rm(a;' — y)) 



L1L2 " 

n 

■f{S\y',d) (D-3) 
where x'^ = ^ and X2 — j^- The additional factor is given by: 

f{x',y',d) ^ exp {lyZl' ~ ^ 

It is independent of n and can be transferred to the right hand side. It gives a 
modified Fourier transform of G{n, lo) for each sector. 

VG(n,c.)exp(27rzn(f' -y*)) = G{x,y,u:)f{x' J ,d) (D-4) 

n 

The next step is to calculate the Fourier transform of G(n, lo) (assuming n to 
be a continuous variable) and use Poissons resummation formula: 



Ztt VfVA V 



(27r)^ ' 2n vpVA 



\x"^y"\ 



KMx"^y"\) 



= Go{x,y,uj) (D-5) 

here: x" = ^ and = ^ for the nodal pair (1,3). For the other pair, vp and 
have to be exchanged. The final expression is obtained by inserting (jP-Sp 
into l(D^ and using jEH): 

G{x,y,u:) = J—Y,G{n,uj)e'^^'"'^''-y'^f{x',y',d) 

n 

= ^f{^+'''^^y^^)Go{x + rn,y,uj) 

= E /V) ( - lUjKoiMx" - y" + m'\) 

+ H^^^^-^S^K,{M{x" - y" + m'\)) (D-6) 
\x" — y"+m'\ / 

The Poisson formula restates the problem in terms of the free space Green's 
functions which are much easier to compute than the ones on a torus. The 
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sums over windings can simply be regularized by omitting the term without 
any winding. This term is divergent and corresponds to the free space Green's 
function. 



E S-wave Contribution 

In several experiments there has been found evidence that the pairing in high- 
Tc Superconductors may not be of pure d-wave but rather a mixture of s- and 
d-wave symmetry. The s-wave contribution prohibits zero-modes: every en- 
ergy eigenvalue is finite, due to the constant energy gap at the Fermi surface. 
Moreover, the s-wave gap introduces further anisotropics between the direction 
perpendicular and the one parallel to the Fermi surface. In the d-wave case, 
the only anisotropy arises from different values for the characteristic velocities 
Vf and va- For the mixed case, the s-wave contribution permits additional 
excitations perpendicular to the Fermi surface but none in the other. Also, the 
energy gap affects the excitations differently. To describe the situation, I replace 
the gap function by: A(/c) —^ A(fc) -I- s(fc), assuming s-wave symmetry for the 
function s{k). It can also be expanded around the nodes to first order: 

sik,+q) = As [{k, + q)l + {k, + q)l^ 

= ^Vskp + Vsqj (E-1) 

For simplicity I neglect vortices, so Vs and va are assumed to be real. After 
inserting (|E-1|) into p^. the Hamiltonian is given by: 

n = nd + -^]ivsidja2 + ^VskFa)-i'j (E-2) 

with a given by: 

■ CTl 

-0-1 

Tid is the pure d-wave Hamiltonian. Although this looks very similar to the 
former case, there are some new aspects to discuss. 

The s-wave contribution does not effect the conservation of the nodon charge. 
Thus, for simplicity I restrict the discussion again to the nodal pair 1/3. The 
most important feature is that the new Hamiltonian docs not have any zero 
modes. There are four distinct energy eigenvalues: 



-El, 2 = ±\J {""fflxY + [vAqy - VsQx + l/^VskF^ 

E^^i = ±^J {vfqxY + {vAQy - VsQx - If-iVskFy (E-3) 

each leading to a unique ground state. It is easy to see, that the energy gap acts 
differently depending on the direction of the excitations. It appears as a mass 
term for the excitations in the x-direction, analogously to the normal s-wave 
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case. However, it only shifts the excitation spectrum in the y-direction. As 
the gap function is constant over the Fermi surface, excitations along the Fermi 
surface see the gap only as an overall shift of energy: a chemical potential. 

I can also be shown that the chiral symmetry is partially broken. The 3- 
direction of the chiral charge is given by the difference of the number of particles 
at node 1 and 3. This conservation is a consequence of the formation of Cooper 
pairs and is not changed by the introduction of an additional s-wavc term. 
However, in the d-wave case, all three components of the chirality were equal 
to each other. I could have chosen any of these to be diagonal. The s-wave 
contribution picks the 3-direction and breaks the symmetry of the other two. 

F Frequency dependence of the Green's func- 
tion 

We want to extract the frequency dependence of the Green's functions after the 
sums over windings are performed: 

ft 

It can be divided into two parts, which I discuss separately: 

(1) = ^ioKoiuV^) (F-1) 

n 

(2) = ^ujKi{ujV^) (F-2) 

n 

For simplicity, I set x = xq, assumed lu > and neglected all constants. Naively, 
you would think that there should not arise any problems as both, (1) and (2) 
are well-defined, finite expressions for each winding. In addition, the Bessel- 
functions decay exponentially for large arguments, thus the sum is truncated 
for large enough winding numbers. In fact, this argument holds for any finite iv. 
However, for small u there are approximately ^ terms for which the argument 
of the Besselfunction has to be assumed as small, that is < 1. So for 

decreasing frequency, the number of terms which which must be summed up 
is growing quadratically in w. Without the alternating factors, this results in 
a -\ behaviour for small oj. The Green's function is ill-defined at w = and 
the contribution of the zero mode to the energy shift has to be calculated in an 
alternative way. In the other sectors, there is no zero mode. So physically, there 
should be no contribution from w = to the energy as there is a low energy 
cut-off. Due to the alternating factors, the terms in the sum can take out each 
other partly and destroy the w — > divergence. 

As the sum cannot be calculated analytically, I approximate it by an integral. 
This is a good approximation, as long as the frequency u) is small and gives the 
correct asymptotic behaviour for a; — > 0. The each alternating factor (—1)" is 
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approximated by cos(n) functions. It is enough to include one of these, as the 
problem becomes spherically symmetric when written as an integral: 



(1) / tPn cos{nx)Ko{iLjn) 



2lT 



d(f) J dnncos{ncos{(j)))KQ(ujn) 
^'^ 1 I cos((/))|arcsinh(cos(0)/aj) 



C0S2 (,/>) + ^2 ^COS2(0)+W 

For the second term, one can derive a similar expression: 



2 



(F-3) 



(2) I d n cos{nx)Ki{ujn) 

r.27r 



d<p / dnncos{ncos{(/)))Ki(ujn) 





Y^C0S2((/)) 

Although the integrals look complicated, they can be solved exactly. Their 
asymptotic behaviour is the same in both cases for small and large to 

(1),(2) ^ ( ^' '^^^^11 (F-5) 
{ ^, uj large 

The approximation we made becomes very bad for large frequencies. In fact, 
for large uj the Green's function shows an exponential decay, instead of the alge- 
braic one in (|F-5p . However, the approximation becomes better the smaller the 
frequencies are. The asymptotic behaviour for small frequencies is, therefore, ac- 
curate. The Green's function approaches a finite constant in the sectors without 
a zero mode. When calculating the energy, one integrates over all frequencies 
and inserts an additional lu in the expression (|76p . In the end, there is indeed 
no contribution to the energy from a; = in absence of a zero mode. The same 
asymptotic behaviour can be expected when the zero mode is excluded from 
the propagator. All remaining modes are gapped and the low-energy behaviour 
should therefore be analogous to the non zero mode sectors. 



G{x,xo,uj) = y2{G{x + n,xo,uj) - d'^y G{y + ft, xo,uj)} (F-6) 

n " 

As the sum is responsible for the w ^ divergence, I neglect the first terms 
and focus on winding numbers where |m| ^ \x\. For large to the difference 
whether you integrate over the torus or whether you choose a specific point 
on the torus becomes negligible and the two terms cancel. By doing a Taylor 
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expansion around to, I show that the subtraetion of the zero mode effectively 
regularizes the Green's function. In fact, it is sufficient to show that Ko{Lo\n\) 
is regularized. The properties for i^i(a;|n| follow out of it. 

The first non-vanishing contribution from the Taylor expansion is second 
order in x. The first derivative, 'VKo{u}\y + n) as well as the mixed term in the 
second order d\d2Ko{w\y + n) vanish because of symmetry reasons. Therefore, 
the first contribution comes from 

{dl + dl)Ko{uj\y + n) 

= {-^K,{u;\n\)+i^^K,{u:\n\)}^^^-^ (F-7) 

Performing the sum, one obtains that the Green's function actually goes to 
zero for a; = which is even more convergent than one would have expected 
from the beginning. This is due to the symmetry on the torus. 
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